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Abstract 


The  focus  of  this  thesis  is  the  experimental  investigation  of  the  dynamics  and  control 
of  a  new  type  of  fast  chaotic  opto-electronic  device:  an  active  interferometer  with 
electronic  bandpass  filtered  delayed  feedback  displaying  chaotic  oscillations  with  a 
fundamental  frequency  as  high  as  100  MHz.  To  stabilize  the  system,  I  introduce 
a  new  form  of  delayed  feedback  control  suitable  for  fast  time-delay  systems.  The 
method  provides  a  new  tool  for  the  fundamental  study  of  fast  dynamical  systems 
as  well  as  for  technological  exploitation  of  chaos. 

The  new  opto-electronic  device  consists  of  a  semiconductor  laser,  a  Mach-Zehnder 
interferometer,  and  an  electronic  feedback  loop.  The  device  offers  a  high  degree  of 
design  flexibility  at  a  much  lower  cost  than  other  known  sources  of  fast  optical 
chaos.  Both  the  nonlinearity  and  the  timescale  of  the  oscillations  are  easily  ma¬ 
nipulated  experimentally.  To  characterize  the  dynamics  of  the  system,  I  observe 
experimentally  its  behavior  in  the  time  and  frequency  domains  as  the  feedback-loop 
gain  is  varied.  The  system  displays  a  route  to  chaos  that  begins  with  a  Hopf  bi¬ 
furcation  from  a  steady  state  to  a  periodic  oscillation  at  the  so-called  fundamental 
frequency.  Further  bifurcations  give  rise  to  a  chaotic  regime  with  a  broad,  flattened 
power  spectrum.  I  develop  a  mathematical  model  of  the  device  that  shows  very 
good  agreement  with  the  observed  dynamics. 

To  control  chaos  in  the  device,  I  introduce  a  new  control  method  suitable  for  fast 
time-delay  systems,  in  particular.  The  method  is  a  modification  of  a  well  known 
control  approach  called  time-delay  autosynchronization  (TDAS)  in  which  the  control 
perturbation  is  formed  by  comparing  the  current  value  of  a  system  variable  to  its 
value  at  a  time  in  the  past  equal  to  the  period  of  the  orbit  to  be  stabilized. The 
current  state  of  a  time-delay  dynamical  system  retains  a  memory  of  the  state  of  the 


system  one  feedback  delay  time  in  the  past.  As  a  result,  the  past  state  of  the  system 
can  be  used  to  predict  the  current  state.  In  order  to  take  advantage  of  this  effect,  the 
new  control  method  forms  a  perturbation  according  to  the  TDAS  scheme  but  delays 
actuation  of  the  control  perturbation  by  a  time  equal  to  the  feedback  delay  time 
of  the  system  to  be  controlled.  This  effectively  sets  the  control-loop  latency  equal 
to  the  feedback  delay  time  of  the  uncontrolled  system.  I  demonstrate  this  control 
method  experimentally  by  stabilizing  a  periodic  orbit  of  the  active  interferometer.  I 
quantify  the  effectiveness  of  the  controller  by  measuring  the  range  of  feedback  loop 
gains  over  which  the  orbit  can  be  stabilized.  The  stabilized  orbit,  which  oscillates 
with  a  frequency  of  51.8  MHz,  is  the  fastest  unstable  periodic  orbit  in  a  chaotic 
system  controlled  experimentally  to  date. 

Application  of  the  new  control  method  requires  the  adjustment  of  two  time 
delays  in  the  controller.  The  first,  the  control  delay  time,  should  equal  the  period 
of  the  orbit  to  be  controlled,  while  the  second,  the  control  loop  latency,  should 
equal  the  feedback  delay  time  of  the  system  to  be  controlled.  I  investigate,  through 
experiments  and  simulations,  the  sensitivity  of  the  method  to  errors  in  setting  these 
time  delays.  I  find  that  the  control  delay  time  must  be  set  exactly  equal  to  the  period 
of  the  orbit  to  minimize  the  control  perturbations  when  the  orbit  is  stabilized.  In 
contrast,  the  control  loop  latency  may  vary  within  a  finite  range  without  affecting 
the  performance  of  the  controller. 
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Chapter  1 

I  ntroduction 

1.1  Chaotic  Dynamics 

Since  efforts  by  the  earliest  astronomers  to  predict  the  motion  of  celestial  bodies, 
dynamics,  or  the  evolution  of  a  system  in  time,  has  been  a  central  concern  of  physical 
science  [1].  Today,  the  study  of  dynamics  is  an  active  branch  of  physics  driven  by 
the  amazing  variety  of  phenomena  observed  in  nonlinear  dynamical  systems,  such 
as  spontaneous  pattern  formation  and  solitons.  Without  question,  one  of  the  most 
widely  studied  nonlinear  dynamical  phenomena  is  chaos.  Chaotic  dynamics  play 
a  role  in  nearly  every  branch  of  physics,  including  condensed  matter  [2,  3],  atomic 
physics  [4],  cosmology  [5,  6],  high  energy  physics  [7,  8],  and  optics  [9,  10]. 

As  the  name  “chaotic”  implies,  these  dynamical  systems  display  highly  irregular, 
even  random-like  oscillations.  In  the  phase  space  of  a  chaotic  system,  trajectories 
beginning  at  initially  close  points  diverge  exponentially  fast.  However,  trajectories 
are  bounded  and  evolve  asymptotically  toward  a  fractal  set  known  as  a  strange  at¬ 
tractor.  Embedded  within  the  strange  attractor  are  a  large  (even  infinite)  set  of 
unstable  periodic  orbits  (UPO).  Loosely  speaking,  these  orbits  form  the  skeleton  of 
the  strange  attractor.  Given  a  knowledge  of  the  UPOs,  one  can  determine  impor¬ 
tant  dynamical  invariants  and  statistical  properties  that  characterize  the  chaotic 
dynamics  [11]- [17]. 

In  1990,  Ott,  Grebogi,  and  Yorke  [18]  demonstrated  that  it  is  possible  to  control 
chaos  in  the  sense  that  unstable  periodic  orbits  in  a  chaotic  attractor  can  be  stabi¬ 
lized  by  feedback  control  applying  only  small  perturbations.  They  pointed  out  that 


1 


since  chaotic  evolution  is  ergodic  a  chaotic  system  will  eventually  come  arbitrarily 
close  to  any  given  UPO.  Thus,  the  controller  simply  waits  for  the  system  to  visit 
the  neighborhood  of  the  desired  orbit  and  then  applies  only  weak  perturbations 
to  keep  it  there.  Numerous  chaos  control  schemes  have  since  been  developed  and 
demonstrated  on  a  wide  variety  of  physical  systems  [19].  Chaos  control  schemes 
are  important  tools  for  building  a  fundamental  understanding  of  chaotic  dynamics 
because  they  provide  a  means  for  identifying  and  characterizing  UPOs. 

Chaos  control  methods  are  also  of  great  practical  importance  given  the  abun¬ 
dance  of  nonlinear  systems  in  nature  and  technology  [20].  These  methods  dramat¬ 
ically  enhance  the  tools  available  to  control  engineers.  Chaos  control  is  now  an 
active  area  of  research  in  the  held  of  control  theory  [21,  22]. 

My  contribution  to  this  held  is  an  experimental  investigation  of  the  control  of 
a  fast  time-delay  chaotic  system.  Delay  dynamical  systems  evolve  in  an  infinite 
dimensional  phase  space  and  display  a  wide  range  of  behavior  including  multista¬ 
bility  and  high- dimensional  chaos.  In  some  respects,  delay-dynamical  systems  are 
analogous  to  spatially  extended  systems  with  the  time  delay  playing  the  role  of 
the  system  size  [23,  24],  While  time-delay  systems  have  been  studied  extensively, 
there  are  many  aspects  of  the  dynamics  of  these  systems  that  are  poorly  under¬ 
stood.  Chaos  control  techniques  suitable  for  time-delay  chaotic  systems  assist  both 
in  the  fundamental  understanding  and  the  technological  exploitation  of  time-delay 
dynamics. 

The  complexity  and  versatility  of  time-delay  dynamics  have  inspired  a  num¬ 
ber  of  proposed  technological  aplications  such  as  secure  communications  [25]- [33], 
dynamic  memory  [34,  35],  and  chaos-based  computation  [36,  37].  Many  of  these 
applications  require  high-speed  chaotic  fluctuations  to  facilitate  rapid  transmission 
or  processing  of  data.  At  the  same  time  they  require  tools  for  manipulating  the 
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chaotic  system  such  as  chaos  control  methods.  However,  control  of  fast  chaotic 
systems  is  particularly  challenging.  The  reason  is  that  all  controllers  require  some 
finite  processing  time,  known  as  control  loop  latency  [38],  between  observing  the 
state  of  the  system  and  perturbing  the  system  to  effect  control.  If  the  state  of  the 
chaotic  system  changes  significantly  during  the  latency  time,  the  control  perturba¬ 
tion  may  not  affect  the  system  in  the  intended  manner  and,  consequently,  control 
may  fail.  Thus,  control  of  fast  chaotic  systems  requires  the  development  of  novel 
control  techniques  that  mitigate  the  influence  of  control  loop  latency. 

1.2  Overview  of  Thesis 

This  dissertation  documents  my  investigation  of  controlling  a  fast  chaotic  time  delay 
system.  It  consists  of  six  chapters.  Chapter  2  presents  a  pedagogical  introduction 
to  control  loop  latency,  and  a  brief  review  of  the  methods  for  controlling  fast  chaos 
developed  prior  to  my  research.  Chapter  3  describes  a  new  time-delay  chaotic  opto¬ 
electronic  device.  The  next  two  chapters  describe  the  application  of  a  new  method 
of  controlling  fast  time-delay  systems  to  the  opto-electronic  device.  The  control 
method  is  a  modification  of  the  time  delay  autosynchronization  (TDAS)  method 
that  is  well-suited  for  fast  time  delay  systems  because  it  allows  for  a  relatively 
large  amount  of  control  loop  latency.  In  Chapter  4,  I  present  an  experimental 
demonstration  of  the  modified  TDAS  control  method.  In  Chapter  5,  I  examine  the 
influence  of  the  control  time  delay  and  latency  on  the  effectiveness  of  control.  I 
summarize  my  main  results  and  discuss  the  direction  of  future  work  in  Ch.  6. 

In  particular,  I  frame  the  problem  of  controlling  fast  dynamics  in  Chapter  2 
through  a  pedogogical  example  and  a  brief  review  of  research  that  has  preceeded 
my  work.  First,  I  examine  the  stabilization  of  the  steady  state  of  a  one  dimensional 
linear  system  for  which  the  domain  of  control  with  and  without  latency  can  be 
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Figure  1.1:  A  new  fast  chaotic  opto-electronic  device:  The  active  interferometer 
with  delayed  bandpass  nonlinear  electronic  feedback. 

determined  analytically.  Then  I  review  several  techniques  for  controlling  chaotic 
dynamics.  Beginning  with  the  original  OGY  method  [18],  I  follow  the  develop¬ 
ment  of  controllers  designed  specifically  for  controlling  fast  chaotic  systems  such  as 
Pyragas’s  time  delay  autosynchronization  [39]. 

In  Chapter  3,  I  introduce  a  flexible  and  inexpensive  new  source  of  fast  optical 
chaos  suitable  both  for  use  in  high-speed  applications  and  low-speed  laboratory 
experiments.  The  system,  shown  schematically  in  Fig.  1.1,  consists  of  a  diode 
laser,  a  Mach-Zehnder  interferometer,  and  an  electronic  delayed  feedback  loop  that 
modulates  the  diode’s  injection  current  according  to  the  power  output  by  the  inter¬ 
ferometer.  I  examine  the  dynamics  of  the  system  when  it  is  running  autonomously 
and  with  external  modulation  provided  by  an  RF  voltage.  In  both  cases,  the  laser 
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emission  undergoes  a  periodic  oscillation  when  the  gain  in  the  feedback  loop  is  low. 
The  frequency  of  oscillation  is  determined  by  the  feedback  delay  time.  As  the  gain 
is  increased,  the  periodic  state  gives  way  to  chaos  through  a  series  of  bifurcations. 
Figure  1.2  shows  a  time  series  of  the  intensity  at  the  interferometer  output  port 


Figure  1.2:  Time  series  of  optical  power  emitted  by  the  interferometer  when  the 
device  is  undergoing  fast  chaotic  oscillation.  Note  that  the  timescale  of  the  oscilla¬ 
tion  is  5-10  ns. 

during  fast  chaotic  oscillation  on  a  nanosecond  timescale.  To  obtain  a  theoreti¬ 
cal  understanding  of  the  behavior  of  the  device,  I  develop  a  detailed  mathematical 
model  that  reproduces  accurately  the  observed  dynamics.  This  fast  time  delay 
system  serves  as  a  testbed  for  the  control  method  introduced  in  the  next  chapter. 

In  Chapter  4,  I  propose  a  new  method  for  controlling  chaos  in  fast  time  delay 
systems.  A  block  diagram  of  such  a  system  is  shown  in  Fig.  4.1.  The  basic  elements 
of  systems  in  this  class  are  a  low-pass  filter,  a  nonlinear  element,  and  a  delayed 
feedback  loop.  Well  known  models  such  as  the  Ikeda  [40]  and  Mackey-Glass  [41] 
equations  are  members  of  this  class  as  well  as  several  chaotic  opto-electronic  devices 
[42]- [45],  [34],  The  new  method  of  control  is  a  modification  of  the  TDAS  control 
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Figure  1.3:  Block  diagram  showing  the  general  form  of  the  class  of  time  delay 
systems  to  which  the  modified  TDAS  method  may  be  applied.  The  basic  elements 
are  a  low  pass  filter,  a  nonlinear  element  (f [x] ) ,  and  a  delayed  feedback  loop.  Well 
known  models  such  as  the  Ikeda  and  Mackey-Glass  equations  are  of  this  form. 

technique  where  the  feedback  is  nonlinear  and  the  control  loop  latency  is  precisely 
set  equal  to  the  time  delay  in  the  uncontrolled  system.  Specifically,  the  form  of  the 
modified  TDAS  feedback  is 

Sp{t)  =  7c  {/  [x  ( t  -  td)\  -  f[x(t-rD-  tr)]}  (1.1) 

where  6p(t)  represents  the  continuous  adjustment  of  an  accessible  system  parameter 
p  about  its  nominal  value,  f  [x  (t)}  is  the  output  of  a  nonlinear  element  in  the  system, 
td  is  the  time  delay  of  the  chaotic  system  to  be  controlled,  7>  is  the  period  of 
the  unstable  orbit  to  be  stabilized,  and  7C  is  the  control  gain.  I  demonstrate  the 
application  of  this  technique  by  using  it  to  control  the  active  interferometer  with 
bandpass  feedback.  To  quantify  the  effectiveness  of  the  controller,  the  domain  of 
control  is  estimated  experimentally  and  numerically. 

In  Chapter  5,  I  continue  my  study  of  the  modified  TDAS  control  method  by 
examining  the  influence  of  the  control  time  delay  and  control  loop  latency  on  the 
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success  of  the  controller.  I  review  studies  of  the  role  of  these  delays  in  the  original 
TDAS  method  and  indicate  aspects  of  these  results  that  should  apply  also  to  the 
modified  TDAS  method.  In  experiments  and  simulations,  the  effectiveness  of  control 
is  observed  as  each  of  these  delays  is  varied  around  its  nominal  value.  The  results 
are  compared  to  expectations  based  on  the  original  TDAS  method. 

In  Chapter  6,  the  final  chapter  of  this  thesis,  I  review  the  main  results  of  my  ex¬ 
periments  and  discuss  the  direction  of  future  work.  The  active  interferometer  with 
bandpass  delayed  feedback  may  be  a  suitable  system  for  a  number  of  possible  ap¬ 
plications  such  as  message  encoding  and  random  number  generation.  The  modified 
TDAS  control  method  represents  a  step  forward  in  the  search  for  effective  control 
methods  for  fast  chaotic  time  delay  systems. 
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Chapter  2 

Control  and  Latency 


Many  techniques  exist  for  “controlling”  dynamical  systems,  that  is,  manipulating 
them  to  produce  a  desired  behavior  [21]- [46].  A  very  common  approach  to  control 
a  system  is  to  use  a  feedback  loop  to  stabilize  an  otherwise  unstable  behavior  of  the 
system  [46].  A  typical  feedback  loop  in  a  controller  consists  of  the  basic  components 
shown  in  Fig.  2.1.  First,  some  information  about  the  state  of  the  system  is  gathered 
through  measurements.  Next,  the  information  is  used  to  determine  how  the  system 
should  be  perturbed  in  order  to  produce  the  desired  behavior  rather  than  the  natural 
evolution  of  the  uncontrolled  system.  Finally,  an  actuator  applies  the  perturbation. 
In  many  cases  these  three  tasks  can  be  completed  in  a  span  of  time  over  which 
the  state  of  the  dynamical  system  changes  very  little.  Then,  for  computational 
convenience,  the  feedback  may  be  considered  to  be  instantaneous.  However,  in  any 
real  feedback  control  scheme,  a  finite  amount  of  time  is  required  for  the  controller 
to  observe  the  system,  determine  the  nature  of  the  perturbation,  and  apply  the 
perturbation. 

This  amount  of  time  is  called  the  control-loop  latency,  labelled  T(  in  Fig.  2.1. 
If  the  unstable  system’s  state  changes  significantly  in  the  time  interval  between 
measurement  and  actuation,  the  perturbation  may  no  longer  be  appropriate  for 
stabilizing  the  system.  Thus,  the  latency  of  a  controller  places  a  practical  limit  on 
the  speed  of  the  instabilities  it  can  stabilize. 

To  get  an  idea  of  when  latency  poses  a  practical  problem,  consider  a  very  simple 
but  commonly  used  control  method  known  as  proportional  feedback  control  [46] .  In 
this  method,  the  controller  produces  a  perturbation  proportional  to  the  difference 


Figure  2.1:  Block  diagram  showing  basic  parts  of  any  feedback  control  scheme. 
The  time,  Tg,  between  measuring  the  state  of  the  system  and  actuating  the  control 
is  known  as  the  control  loop  latency. 

between  a  measurement  of  the  state  of  the  system  and  a  reference  value.  Suppose 
the  state  of  the  system  is  quantified  by  a  voltage  (produced  by  a  measurement 
tranducer,  for  example).  First,  consider  a  controller  implemented  using  a  personal 
computer  with  an  analog-to-digital  converter  to  record  the  measurement,  some  high- 
level  software  like  LabView  to  evaluate  the  difference  between  the  measured  state 
and  the  reference  value,  and  a  digital-to-analog  converter  to  actuate  the  control.  In 
that  case,  the  latency  would  be  on  the  order  of  10  /is.  Therefore,  unstable  dynamics 
involving  frequencies  beyond  about  50  kHz  would  be  too  fast.  If  the  controller  is 
implemented  using  analog  electronics  such  as  operational  amplifiers  (10  MHz  typical 
bandwidth  at  unity  gain),  then  instabilities  with  a  ~  5  MHz  characteristic  frequency 
could  likely  be  stabilized.  Using  RF  electronics,  which  typically  have  a  bandwidth 
on  the  order  of  1  GHz,  latency  could  be  reduced  to  a  few  nanoseconds,  allowing 
control  of  ~  500  MHz  fluctuations. 
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In  the  case  of  proportional  feedback  control,  any  one  of  these  implementations 
is  possible.  However,  such  a  simple  controller  is  not  effective  in  many  cases,  such 
as  stabilizing  unstable  periodic  orbits  in  chaotic  or  other  nonlinear  systems.  More 
complex  controllers  that  produce  very  good  results  when  applied  to  slow  systems 
cannot  be  implemented  so  easily  at  higher  frequencies.  Thus,  there  is  a  clear  need 
for  control  techniques  designed  specifically  to  work  in  the  regime  of  “fast”  dynamics, 
i.e.  when  the  latency  is  comparable  to  the  characteristic  time-scale  of  the  dynamical 
system.  The  following  chapters  will  describe  a  new  approach  to  controlling  one 
class  of  fast  chaotic  systems,  but  first  a  review  of  the  basic  problem  of  control  with 
latency  and  past  approaches  to  solving  it  is  presented  here.  In  Sec.  2.1,  I  illustrate 
the  problem  posed  by  control  loop  latency  through  a  simple  example  for  which  the 
domain  of  control  and  it’s  dependence  on  latency  can  be  calculated  analytically.  In 
Sec.  2.2,  I  review  the  development  of  controllers  for  fast  chaotic  systems. 

2.1  Control-Loop  Latency:  A  Linear  Example 

In  this  section,  I  introduce  several  ideas  needed  to  understand  the  problem  of  latency 
in  a  controller.  The  concepts  of  stability,  control,  and  latency  are  illustrated  by  a 
the  example  of  a  first-order  linear  dynamical  system.  An  analysis  of  this  system 
provides  a  simple  demonstration  of  how  control  fails  due  to  latency. 

Consider  the  very  simple  dynamical  system  described  by  the  linear  differential 
equation 

x  =  ax.  (2.1) 

The  point  x  =  0  is  the  only  fixed  point  of  the  system,  that  is,  a  point  in  phase 
space  where  x  —  (1  [48].  Suppose  the  fixed  point,  denoted  by  x*,  is  the  desired 
behavior  of  the  system.  Consider  the  effect  of  a  small  perturbation  from  the  fixed 
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point  <5  X  X%  .  The  evolution  of  this  perturbation  is  obtained  by  substituting 
into  Eq.  2.1  to  obtain  the  differential  equation 

6  =  aS,  (2.2) 

the  solution  to  which  is 

S(t)  =  SQeat,  (2.3) 

where  6a  is  the  initial  perturbation  to  the  system.  If  a  <  0,  the  perturbation  decays 
and  the  fixed  point  is  stable.  On  the  other  hand,  if  a  >  0,  the  perturbation  will 
grow  so  the  fixed  point  is  unstable.  In  a  situation  where  the  fixed  point  represents  a 
desired  behavior  but  is  unstable,  a  controller  can  be  added  to  the  system  to  maintain 
the  desired  behavior.  I  will  assume  a  >  0  for  the  rest  of  this  section. 

Note  that  a  is  the  exponential  growth  rate  of  perturbations  and  therefore  is  the 
natural  time  scale  of  the  dynamical  system.  When  control  is  added  to  the  system  it 
is  reasonable  to  expect  that  latency  will  degrade  the  performance  of  the  controller 
when  it  is  of  the  order  of  or  greater  than  a-1.  In  this  case,  the  perturbation  can  grow 
significantly  before  the  controller  can  respond.  Before  demonstrating  this  formally, 
I  will  examine  the  case  of  instantaneous  feedback. 

To  stabilize  the  fixed  point,  I  apply  proportional  feedback  control  as  described 
above.  The  system  state  is  compared  to  a  reference  state  (the  fixed  point  in  this 
case)  and  a  signal  proportional  to  the  difference  is  fed  back  to  the  system  affecting 
its  future  evolution.  When  the  feedback  is  instantaneous,  the  dynamical  system 
plus  controller  is  described  by  the  differential  equation 

x  =  ax  +  b(x  —  x*),  (2.4) 

where  b  denotes  the  strength  of  the  feedback  or  control  gain.  Since  x*  =  0,  Eq.  2.4 
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can  be  rewritten  as 


x  =  (a  +  b)x.  (2.5) 

To  determine  the  controller’s  effect  on  the  stability  of  the  fixed  point,  I  again  con¬ 
sider  a  small  perturbation,  6  =  x  —  x*.  Equation  2.5  implies  the  perturbation  will 
evolve  according  to  the  differential  equation 

8  =  (a  +  b)8  (2.6) 

whose  solution  is  given  by 

8(t)  =  80e^a+b)t.  (2.7) 

From  Eq.  2.7,  the  effect  of  the  controller  is  clear.  If  b  <  —a,  the  perturbation 
decays  and  the  fixed  point  is  stable.  A  useful  method  of  visualizing  the  effect  of  a 
controller  is  to  plot  the  control  gain  b  versus  a  system  parameter,  such  as  a.  The 
region  in  this  plot  where  the  unstable  state  is  stabilized  by  the  controller  is  referred 
to  as  the  domain  of  control.  Figure  2.2  shows  the  domain  of  control  for  the  example 
system  with  proportional  feedback  control  and  no  control-loop  latency. 

Having  demonstrated  how  a  desired  behavior  of  a  dynamical  system  is  attained 
through  feedback  control,  two  important  points  should  be  noted.  First,  in  this 
example,  the  desired  behavior  happened  to  be  a  fixed  point  or  steady  state.  In 
general,  unstable  periodic  states  (more  accurately  unstable  periodic  orbits  or  UPOs) 
in  phase  space  can  be  controlled,  too  [21,  22],  In  fact,  research  on  controlling  chaos 
tends  to  focus  on  stabilizing  UPOs  rather  than  fixed  points.  The  reason  for  this 
is  that  chaotic  systems  typically  contain  an  infinite  number  of  UPOs,  presenting  a 
surprising  degree  of  flexibilty  from  a  single  system.  Second,  adding  the  controller 
to  the  dynamcal  system  does  not  change  the  location  of  the  fixed  point;  only  the 
stability  of  the  fixed  point  is  altered.  Some  schemes  for  control  do  involve  making 
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Figure  2.2:  Domain  of  control  of  the  fixed  point  at  origin  of  the  first-order  linear 
dynamical  system  with  instantaneous  proportional  feedback  control.  The  axes  show 
the  bifurcation  parameter,  a,  and  the  control  gain,  b. 

dramatic  changes  to  the  phase  space  of  the  unstable  system  such  as  altering  UPOs  or 
creating  new  UPOs.  However,  it  is  often  desirable  that  the  behavior  of  the  unstable 
system  should  be  preserved,  e.  g.  to  minimize  the  power  used  by  the  controller. 
Thus,  the  controller  should  not  introduce  new  UPOs  into  the  dynamics. 

To  investigate  the  effect  of  latency  on  the  controller,  I  determine  how  the  domain 
of  control  changes  as  latency  becomes  significant.  Suppose  the  control  signal  is 
applied  at  a  time  T£  later  than  it  would  be  if  it  were  instantaneous.  The  evolution 
of  the  dynamical  system  is  then  given  by 

x(t)  =  ax(t)  +  b{x{t  —  Ti)  —  a;*}, 

=  ax(t)  +  bx(t  —  Tg).  (2.8) 

Once  again,  the  controller’s  effect  on  the  stability  of  the  fixed  point  is  determined 
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by  examining  the  evolution  of  a  small  perturbation  <5  =  x  —  x*.  Equation  2.8  implies 
the  perturbation  will  evolve  according  to 

6(t)  —  a6(t)  +  b6(t  —  Tg)  (2.9) 

To  solve  this  equation,  I  insert  an  exponential  trial  solution  into  Eq.  2.9  of  the  form 

6  (t)  =  60ext  (2.10) 

where  A  is  the  eigenvalue  of  the  fixed  point.  The  result  is  the  characteristic  equation 

A  =  a  +  be~Xr  , 

A  -a-be-Xr  =  0.  (2.11) 

The  left  hand  side  of  Eq.  2.11  is  known  as  the  characteristic  quasipolynomial  and 
will  be  denoted  $(A).  If  the  fixed  point  is  to  be  stable,  i.e.  the  perturbation  is 
to  decay,  there  must  be  no  solutions  to  Eq.  2.11  containing  a  positive  real  part. 
Equivalently,  <b(A)  must  have  no  roots  with  positive  real  parts.  The  region  in  the 
ab  parameter  space  where  this  condition  is  satisfied  is  the  domain  of  control. 

To  End  this  region,  I  use  the  method  of  D-partition  [47].  I  divide  the  ab  plane 
into  distinct  regions  separated  by  curves  on  which  <3>(z)  has  at  least  one  root  with 
real  part  equal  to  zero.  At  all  points  within  one  such  region  of  the  plane,  $(A)  has 
the  same  number  of  roots  with  a  positive  real  part.  This  partition  of  the  parameter 
space  is  known  as  a  D-partition.  To  locate  the  domain  of  control,  I  identify  the 
particular  region  in  which  that  number  is  zero.  The  first  boundary  curve  is  found 
by  setting  A  equal  to  zero  to  get  the  line 

a  =  -b.  (2.12) 

To  obtain  the  rest,  assume  A  =  iy  so  Eq.  2.11  becomes 

iy-a-be~ivr  =0.  (2.13) 
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Figure  2.3:  H-partition  of  the  ab  parameter  space.  The  dashed  lines  show  curves 
on  which  the  characteristic  quasipolynomial  has  at  least  one  root  with  a  real  part 
equal  to  zero. 


Separating  the  real  and  imaginary  parts  gives  the  parametric  form  of  the  infinite 
set  of  curves  that  make  up  the  remaining  boundaries 


cos  rey 

a  =  y- - , 

sm  Tty 


(2.14) 


b= - — .  (2.15) 

sm  rey 

One  curve  defined  by  these  equations  meets  the  line  a  =  —b  at  a  cusp  point  (l/iy, 
—  1/t^).  The  first  few  bounding  curves  closest  to  the  origin  are  shown  in  Fig.  2.3. 
A  closer  view  of  the  region  surrounding  the  origin  is  given  in  Fig.  2.4 

In  the  simple  case  where  a  <  0  and  no  control  is  applied  (6  =  0),  the  fixed  point 
is  stable  so  no  solutions  to  Eq.  2.11  have  positive  real  parts.  Therefore,  for  all 
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Figure  2.4:  .D-partition  of  the  ab  parameter  space  used  to  locate  the  domain  of 
control.  In  Region  I,  no  roots  of  the  characteristic  quasipolynomial  have  positive  real 
parts.  In  Regions  II  and  III,  one  and  two  roots  have  positive  real  parts,  respectively. 
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points  in  Region  I  of  Fig.  2.4  there  are  no  solutions  to  Eq.  2.11  with  positive  real 
parts.  To  determine  the  number  of  roots  with  positive  real  part  in  Region  II,  note 
that  along  the  line  6  =  0  the  characteristic  equation  reduces  to  the  simple  form 

A  =  a,  (2.16) 

so  there  is  clearly  only  one  root  with  a  positive  real  part.  (As  b  — >  0,  the  real  parts 
of  all  the  roots  except  one  approach  —  oo.) 

To  determine  the  number  of  roots  with  positive  real  part  in  region  III,  consider 
the  sign  of  the  differential  of  the  root  of  <h(A)  with  zero  real  part  as  a  boundary  is 
crossed.  Specifically,  if 

$  (A,  a,  b)  —  0,  (2.17) 


then 


<9d>  <9d>  n 

— — aA  H — — — da  H — —do  —  0. 
oX  da  db 


(2.18) 


The  differential  of  the  real  part  of  the  root  is 

Re  (dX)  =  Re 


—  H da  —  | ^db 
d\ 


(2.19) 


Moving  from  Region  II  to  III  across  the  line  a  =  —b,  assuming  da  <  0 ,  db  =  0,  and 


b  <  —  1/t^,  Eq.  2.19  becomes 

Re  (dX)  =  d° l  >  0.  (2.20) 

V  7  1  +  bn  V  7 

The  real  part  receives  a  positive  increment  implying  that  points  in  Region  III  have 
at  least  one  more  root  with  positive  real  part  than  in  Region  II.  A  similar  analysis 
of  the  other  boundaries  shows  that  further  roots  with  positive  real  parts  appear  as 
each  boundary  is  crossed  on  the  line  a  =  0  moving  away  from  Region  I.  Therefore, 
Region  I  is  the  only  region  with  no  unstable  roots. 
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The  domain  of  control  is  simply  the  section  of  Region  I  to  the  right  of  the  line 
a  =  0  as  shown  in  Fig.  2.5.  In  contrast  to  the  latency-free  case,  no  control  is 


Figure  2.5:  Domain  of  control  of  the  first-order  linear  system  with  proportional 
feedback  control  and  non-zero  control  loop  latency.  Control  is  only  possible  when 
latency,  ti,  is  less  than  the  response  time  of  the  system,  a _1. 


possible  when  T£  <  1/a.  Control  can  only  be  achieved  if  the  latency  is  shorter  than 
the  characteristic  timescale  of  the  system.  This  result  is  consistent  with  the  intuitive 
argument  presented  above  that  perturbations  to  the  system  may  grow  too  rapidly 
for  the  control  response  to  be  effective.  Figure  2.5  also  shows  that  the  domain  of 
control  is  of  finite  extent  even  when  <  1/a,  whereas  it  extends  indefinitely  when 
latency  is  not  present. 

Analogous  effects  of  control  loop  latency  have  also  been  found  in  chaos  control 
schemes.  Sukow  et  al.  [49]  investigated  the  effect  of  latency  on  control  of  a  fast 
chaotic  electronic  circuit  using  two  different  control  schemes  known  as  time-delay 
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auto  synchronization  or  TDAS  and  extended  time-delay  autosynchronization  or  ET- 
DAS  (an  explanation  of  these  schemes  appears  in  the  next  section).  In  both  cases, 
the  domain  of  control  decreased  in  size  as  the  latency  was  increased  until  control 
was  finally  lost.  The  maximum  latency  at  which  control  was  attained  was  as  much 
as  four  times  larger  with  ETDAS  than  TDAS,  but  ETDAS  control  failed  when  the 
latency  reached  ~  86%  of  the  correlation  time  of  the  uncontrolled  orbit  even  in  the 
best  case  observed.  Just  et  al.  [50]  developed  an  approximate  prediction  for  the 
critical  latency  at  which  TDAS  control  fails.  They  predict  TDAS  control  can  be 
achieved  when 

(l  _  "zeA 

Tie  <  (2.21) 

UTp 

where  t>  and  v  are  the  period  and  Floquet  exponent  (or  average  growth  rate  of 
perturbations)  of  a  the  UPO  to  be  stabilized,  respectively.  This  prediction  was 
tested  in  experiments  on  a  nonlinear  electronic  circuit  with  Rossler  type  behavior. 
Control  failed  experimentally  when  the  latency  reached  a  value  ~  11%  of  period  of 
the  UPO.  Equation  2.21  predicted  failure  at  ~  12.5%  of  the  period,  in  reasonable 
agreement  with  the  experiemental  results. 

2.2  Latency  and  Controlling  Fast  Chaos 

The  dynamical  system  examined  in  the  previous  section  shows  that  control  of  even 
the  simplest  fixed  point  may  fail  in  the  presence  of  control-loop  latency.  Control 
engineers  have  contended  with  the  problem  of  latency  for  decades  [51] .  Most  work  on 
this  problem  has  focused  on  the  properties  fixed  points  of  linear  dynamical  systems. 
For  example,  methods  for  determining  stability,  controllability  and  robustness  to 
latency  of  linear  systems  have  been  developed  (see  Refs.  [52,  53]  and  the  references 
therein).  In  recent  years,  physicists  have  taken  note  of  this  problem  as  it  presents  a 
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practical  limit  on  the  application  of  chaos-control  schemes  to  stabilize  UPOs  of  fast 
chaotic  systems.  High-speed  chaotic  optical  systems  have  been  studied  that  have 
characteristic  timescales  ~  1  ns  or  less  [54]- [57].  The  possibility  of  employing  fast 
controlled  or  synchronized  chaotic  optical  systems  in  novel  communication  systems 
[25]-  [33]  has  enhanced  the  urgency  of  finding  new  techniques  for  dealing  with 
latency.  In  this  section,  I  will  trace  the  evolution  of  fast  chaos-control  techniques 
towards  controllers  that  require  less  computation  in  order  to  minimize  latency. 

The  first  technique  for  controlling  chaotic  systems  was  proposed  by  Ott,  Grebogi, 
and  Yorke  [18]  in  1990  and  is  commonly  known  as  the  OGY  method.  According 
to  this  method,  the  dynamical  system  to  be  controlled  is  represented  by  a  discrete 
map  of  the  form 

Zi+l  =  F  (Zj ,  p) ,  (2.22) 

where  Zj  is  a  vector  representing  the  system  state  ( e.g .,  the  current  location  on  a 
surface  of  section  in  a  time-delay  reconstructed  phase  space)  and  p  is  an  accessible 
system  parameter  that  can  be  varied  about  a  nominal  value  p.  The  UPO  to  be 
stabilized  is  an  unstable  fixed  point,  denoted  Z*,  of  the  uncontrolled  system,  i.e. 
the  map  F  (Z|,  p).  Time  series  analysis  is  used  to  determine  a  linear  approximation 
of  the  dynamics  near  the  fixed  point  of  the  form 

Zi+i  -  Z*  =  A  (Zj  -  Z*)  +  B  (p  -  p),  (2.23) 

where  A  represents  the  system  response  to  deviations  from  the  fixed  point  of  the 
uncontrolled  system  and  B  represents  the  response  to  deviations  in  the  parameter, 
p,  respectively.  Then  feedback  is  implemented  by  varying  the  parameter  p  around 
a  nominal  value  p  according  to  the  rule 

P-P  =  -KT  [Zj  -Z*],  (2.24) 
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where  the  matrix  K  determines  the  strength  of  the  feedback.  Substituting  Eq.  2.24 
into  Eq.  2.23,  I  obtain 

Zj+1  -z*  =  (A  -  BKT)(Zi  -z*).  (2.25) 

The  fixed  point  is  stable  if  the  eigenvalues  of  the  matrix  A  —  B  K T  have  modu¬ 
lus  smaller  than  unity.  A  well  known  technique  from  control  engineering  called 
pole  placement  [46]  is  used  to  determine  the  appropriate  value  of  K  to  ensure  this 
condition  is  satisfied  [59]. 

The  effectiveness  of  the  OGY  method  was  demonstrated  experimentally  by  Ditto 
et  al.  [60]  who  applied  it  successfully  to  a  chaotic  magnetoelastic  ribbon  in  an  ap¬ 
plied  AC  magnetic  field.  A  personal  computer  was  used  to  perform  time  series 
analysis  to  obtain  the  necessary  linear  approximations  and  vector  calculations  in- 
vloved  in  evaluating  Eq.  2.24.  As  stated  earlier,  this  introduced  many  microseconds 
of  latency  into  the  control  loop.  This  was  not  a  problem  in  the  work  of  Ditto  et 
al.  because  the  characteristic  timescale  of  the  oscillations  they  suppressed  was  0.85 
seconds  (the  period  of  the  applied  AC  field).  However,  such  latency  is  a  problem 
when  the  system  to  be  controlled  oscillates  at  frequencies  above  ~  10  kHz. 

An  early  modification  of  the  OGY  method  was  occasional  proportional  feedback 
[61].  In  this  scheme,  a  scalar  state  variable  z(t)  of  the  dynamical  system  to  be 
controlled  is  sampled  at  regular  intervals  to  create  a  discrete  representation  of  the 
system.  If,  on  the  ith  sample,  the  observed  value  of  the  state  variable  z,  falls  within 
a  window  of  width  W  centered  on  the  value  z*,  control  is  activated  by  adjusting  an 
accessible  parameter  p  from  it’s  nominal  value  p  by  an  amount 

Spi  =  7  (zi-  z*) ,  (2.26) 

where  7  is  the  control  gain.  Otherwise,  the  parameter  p  is  maintained  at  it’s  nominal 
value  p.  The  appropriate  values  of  the  control  gain  7,  window  size  W,  and  center 
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location  z*  can  quickly  be  obtained  by  trial  and  error.  The  computationally  costly 
linear  fitting  and  vector  calculations  involved  in  the  OGY  method  are  avoided  in 
OPF.  As  a  result,  this  technique  is  easily  implemented  using  analog  electronics  that 
produce  a  smaller  latency.  Roy  et  al.  [62]  employed  the  fastest  implementation  of 
OPF  reported  in  the  literature  so  far  to  stabilze  chaotic  oscillations  of  a  multimode 
laser  with  characteristic  frequency  118  kHz. 

Myneni  et  al.  [63]  simplified  the  OPF  technique  to  further  reduce  control-loop 
latency.  In  their  scheme,  the  system  state  z(t)  (scalar  or  vector)  is  observed  contiu- 
ously.  When  the  system  enters  a  predefined  window  W  in  phase  space,  an  accesible 
parameter  p  is  perturbed  from  it’s  nominal  value  p  by  an  amount  6p  of  fixed  mag¬ 
nitude  and  sign.  The  control  perturbation  is  turned  off  when  the  system  exits  from 
W.  The  boundary  and  location  of  W  are  chosen  so  that  the  intersection  of  W 
and  the  desired  UPO  is  zero  and  the  transit  time  through  W  increases  along  the 
local  unstable  directions.  As  in  the  OPF  method,  the  apropriate  window  size  and 
location  and  the  magnitude  and  sign  of  6p  can  be  determined  by  trial  and  error. 
Myneni  et  al.  implemented  this  scheme  using  high-speed  electronic  components  on 
a  custom  printed  circuit  board  achieving  a  control-loop  latency  of  just  4.4  ns.  With 
this  controller,  they  successfully  stabilized  UPOs  of  a  chaotic  Colpitts  oscillator 
with  a  characteristic  frequency  of  19.1  MHz. 

The  techniques  described  so  far  represent  one  line  in  the  evolution  of  fast  chaos 
controllers.  A  distinct  line  of  control  techniques  is  based  on  continuous  feedback 
proportional  to  a  past  state  of  the  system.  Unlike  the  controllers  mentioned  previ¬ 
ously,  these  techniques  require  no  switching  of  the  perturbation  from  one  value  to 
the  next  or  on  and  off.  The  process  of  switching  must  be  done  on  a  much  shorter 
time  scale  than  the  chaotic  oscillations.  Therefore,  in  principle  this  approach  can 
be  applied  to  systems  that  are  so  fast  that  suitable  switching  technology  is  not  ex- 
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ist  such  as  high-speed  optical  systems.  In  theory,  at  least,  many  of  these  delayed 
feedback  techniques  can  even  be  implemented  wholly  in  the  optical  domain. 

The  first  such  method  was  time-delay  autosynchronization  or  TDAS  [39].  In 
this  method,  the  dynamical  system  of  interest  is  subjected  to  continuous  feedback 
through  an  accessible  parameter  p  of  the  form 

=  ~l{z{t)  -z{t-TP)},  (2.27) 

where  z(t)  is  an  observed  system  variable,  7  is  the  control  gain,  and  r  is  fixed  a 
time  interval.  Control  is  achieved  by  setting  r  equal  to  the  period  of  the  UPO  to 
be  stabilized  and  choosing  an  appropriate  control  gain.  If  no  accurate  model  of 
the  unstable  system  is  available,  the  control  gain  and  delay  can  simply  be  swept 
to  locate  suitable  values.  However,  for  situations  where  a  model  is  available,  there 
exist  numerical  [64]  and  analytical  methods  [65,  66]  for  predicting  the  domain  of 
control. 

Time-delay  autosynchronization  has  been  successfully  applied  to  such  diverse 
experimental  systems  as  electronic  circuits  [67,  49],  Taylor-Couette  fluid  flow  [68], 
an  15NH3  laser  ,  and  plasma  instabilities  [70,  71].  Sukow  et  al.  [49]  reported  the 
fastest  implementation  of  TDAS  prior  to  the  current  study  using  it  to  stabilize  a  10.1 
MHz  diode  resonator.  Their  controller,  constructed  using  high-frequency  analog 
electronics,  had  a  latency  of  10  ns  [38].  However,  due  to  the  computational  simplicity 
of  TDAS,  much  faster  implementations  certainly  seem  possible.  For  example,  an 
all-optical  implementation  of  TDAS  has  been  proposed  and  successfully  applied  to 
a  numerical  model  of  a  two- level  laser  system  [73]. 

Some  interesting  generalizations  of  TDAS  have  been  investigated.  Just  et  al.  [66] 
considered  the  possibility  of  using  a  nonlinear  function  g  [z{t)\  of  a  system  variable 
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to  create  the  feedback  signal 


8p(t)  =  -7  {g[z{t)]  ~  g[z(t  -  tp)]} 


(2.28) 


in  experimental  situations  where  the  variable  z(t)  itself  is  not  directly  accessible. 
Nakajima  and  Ueda  showed  that  for  some  orbits  where  TDAS  fails  control  can  be 
achieved  by  setting  the  control  time  delay  to  half  the  period  of  the  orbit  [74] .  Socolar 
et  al.  [72]  introduced  a  generalization  of  TDAS  where  more  than  one  past  state  of 
the  system  is  used  to  create  the  feedback  signal.  In  this  scheme,  called  extended 
time-delay  autosynchronization  or  ETDAS,  the  control  signal  is  of  the  form 


6p(t) 


-7 


z(t)  —  (1  —  R)  Rk  1z(t  —  krp ) 

k=  1 


(z(t)  -  z(t  -  Tp)) 


(2.29) 


where  0  <  R  <  1  determines  the  weight  of  past  states  in  the  sum.  This  technique 
has  been  demonstrated  experimentally  on  electronic  circuits. 

In  Ch.  4,  I  introduce  a  new  modification  of  TDAS  suitable  to  fast  chaotic 
systems  with  a  time  delay.  Many  such  systems  are  currently  being  studied  for 
use  in  high  data  rate  chaos-based  communications  systems.  The  following  chapter 
describes  one  such  system  in  detail.  Applying  control  to  fast  time-delay  systems 
will  almost  certainly  mean  dealing  with  the  issue  of  latency.  However,  the  current 
state  of  a  time-delay  dynamical  system  retains  a  memory  of  the  state  of  the  system 
one  feedback  delay  time  in  the  past.  As  a  result,  the  past  state  of  the  system  can 
be  used  to  predict  the  current  state.  The  modified  TDAS  control  scheme  takes 
advantage  of  this  effect  by  delaying  actuation  of  the  control  perturbation  by  a  time 
equal  to  the  feedback  delay  time  of  the  system  to  be  controlled. 
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Chapter  3 

The  Active  I  nterferometer  with 
Bandpass-Filtered  Time-Delayed 
Feedback 


Dynamical  systems  with  a  time-delay  or  memory  effect  can  display  an  astonishing 
range  of  complex  behavior.  For  example,  Ikeda  et  al.  [75]  studied  a  passive  optical 
resonator  described  by  the  equation 

x(t)  =  —x(t)  +  n  sin  [x(t  —  tit)]  ,  (3.1) 

where  //  and  tn  are  the  strength  of  nonlinearity  and  the  time-delay,  respectively. 
With  no  time-delay  (£r  =  0)  this  one  dimensional  dynamical  system  is  limited  to 
simple  fixed-point  dynamics.  However,  when  tn  »  1  this  system  can  exhibit  multi- 
stability  and  high- dimensional  chaos  [40].  Such  complex  time-delay  dynamics  arise 
in  many  physical,  biological,  chemical,  and  engineering  models  [52],  [76]- [78].  Opti¬ 
cal  time-delay  systems,  in  particular,  have  been  a  subject  of  continued  interest  over 
the  past  two  decades.  Besides  the  passive  optical  resonator  [40,  76]  mentioned  above, 
other  important  systems  include  semiconductor  lasers  with  optical  or  electro-optical 
feedback  [79]-  [81],  CO2  lasers  with  delayed  feedback  via  an  intracavity  electro-optic 
modulator  [82,  83],  and  electro-optic  hybrid  devices  with  delayed  feedback  [34], 
More  recently,  the  possibility  of  chaos-based  technology  has  sparked  interest  in 
fast  chaotic  optical  time-delay  systems.  Proposed  technological  applications  of  such 
systems  include  dynamic  memory  [34],  dynamics-based  computation  [36,  37],  and 
communications  [25]- [33].  In  each  of  these  applications,  high-speed  chaotic  fluc¬ 
tuations  are  needed  to  facilitate  rapid  processing  or  transfer  of  information.  For 
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example,  chaos-based  secure  communication  at  250  Mbit/s  has  been  achieved  us¬ 
ing  an  erbium-doped  fiber  ring  laser  (EDFL)  as  a  source  of  picosecond  chaotic 
oscillations  [29].  The  complexity  of  time-delay  dynamical  systems  may  also  be  an 
advantage  in  these  applications.  Research  on  the  security  of  chaos-based  commu¬ 
nication  schemes  indicates  high- dimensional  dynamics  are  necessary  to  ensure  that 
nonlinear  forecasting  techniques  cannot  be  used  by  an  eavesdropper  to  extract  a 
message  from  the  chaotic  carrier  [88]-  [90]. 

In  this  chapter,  I  introduce  a  new  fast  chaotic  time-delay  dynamical  system  as  a 
flexible,  inexpensive  source  of  fast  optical  chaos.  Both  the  timescale  and  the  com¬ 
plexity  of  the  dynamics  are  easily  tuneable.  Compared  to  other  time-delay  optical 
system  such  as  the  EDFL  mentioned  above,  this  system  is  constructed  with  inex¬ 
pensive,  readily  accessible  components.  These  characteristics  make  it  an  attractive 
experimental  system  for  studying  time-delay  dynamics  and  control.  In  the  next 
chapter,  the  device  will  serve  as  a  testbed  for  a  new  chaos  control  method. 

Figure  3.1  shows  a  schematic  diagram  of  the  device.  Light  emitted  from  a 
laser  interferometer  falls  on  a  photodiode.  The  electrical  signal  from  the  diode  is 
amplified,  propagates  through  a  long  coaxial  cable,  and  is  used  to  modulate  the 
injection  current  of  the  laser.  The  propagation  time  through  the  cable  is  long 
compared  to  the  response  time  of  the  laser  interferometer  itself.  The  nonlinearity  of 
the  interferometer  coupled  with  the  delay  in  the  feedback  loop  combine  to  produce 
a  range  of  steady  state,  periodic,  and  chaotic  behavior. 

This  new  device  has  some  important  differences  with  previous  optical  time-delay 
systems.  First,  Ohtsubo  et  al.  [85]-  [43]  studied  an  active  interferometer  with  elec¬ 
tronic  feedback,  but  used  low  frequency  limiting  the  oscillation  frequency  to  tens 
of  kilohertz.  In  my  system,  the  high-frequency  electronic  feedback  signal  is  ap¬ 
plied  to  a  semiconductor  laser  through  a  radio- frequency  device  called  a  bias-T 
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F  igure  3.1:  Schematic  of  active  interferometer  with  delayed  band-pass  filtered  elec¬ 
tronic  feedback.  The  device  produces  chaotic  oscillations  on  nanosecond  timescales. 

allowing  me  to  achieve  >100  MHz  oscillations.  Using  state  of  the  art  technologies, 
10-100  GHz  frequencies  should  be  attainable.  Second,  the  low-frequency  electron¬ 
ics  in  Ohtsubo’s  system  produced  a  low  pass  characterisitic  in  the  feedback  loop. 
In  contrast,  the  bias-T  used  in  my  system  produces  a  bandpass  characteristic  in 
the  feedback  loop  blocking  any  DC  component  from  circulating.  The  dynamics  of 
time-delay  systems  with  band-pass  filtered  feedback  is  a  largely  unexplored  topic. 
However,  it  is  of  increasing  importance  for  fast  chaos  applications  since  many  high 
speed  opto-electronic  devices  are  AC-coupled.  Two  interesting  results  that  have 
been  reported  are  the  use  of  bandpass  filtering  to  tailor  the  spectrum  of  the  chaotic 
signal  to  fit  a  desired  communication  band  and  the  observation  of  apparently  very 
high-dimensional  chaos  [44,  87]. 
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In  Sec.  3.1,  I  describe  the  experimental  system.  In  Sec.  3.2,  I  develop  a  theoret¬ 
ical  model  of  the  active  laser  interferometer  with  bandpass  feedback.  In  Sec.  3.3,  I 
describe  the  measurements  necessary  to  determine  the  parameters  of  the  model.  In 
Sec.  3.4,  the  dynamics  of  this  new  system  are  explored  through  a  survey  of  the  route 
to  chaos.  Section  3.5  describes  this  effect  of  external  modulation  on  the  system.  In 
Sec.  5.3  I  discuss  the  issues  involved  in  scaling  the  system  to  run  at  faster  speeds, 
a  necessity  for  use  in  a  real  communication  system. 

3.1  Experimental  Apparatus 

In  this  section  I  describe  my  experimental  implementation  of  the  active  interferom¬ 
eter  with  bandpass-filtered  delayed  feedback.  The  light  source  is  a  0.65  [irri  wave¬ 
length  AlGalnP  diode  laser  (Hitachi  HL6501MG)  with  a  multi-quantum  well  struc¬ 
ture.  The  diode  is  housed  in  a  commercial  mount  (Thorlabs  TCLDM9)  equipped 
with  a  bias-T  for  adding  an  RF  component  to  the  injection  current,  and  tem¬ 
perature  regulation  circuitry.  Thermoelectric  coolers  connected  to  a  proportional- 
integral-derivative  feedback  controller  (Thorlabs  TEC2000)  provide  1  rnK  temper¬ 
ature  stability,  minimizing  frequency  and  power  drift  due  to  heating  effects.  The 
light  generated  by  the  laser  is  collimated  by  a  (Thorlabs  C230TM-B,  f=4.5mm) 
lens,  producing  an  elliptical  beam  (1  mm  x  5  mm)  with  a  maximum  output  power 
of  35  niW  at  the  maximum  operating  current  of  90  mA. 

The  laser  beam  is  directed  into  a  Mach-Zehnder  interferometer.  The  interferom¬ 
eter  contains  two  paths  for  the  light  to  follow  whose  lengths  differ  by  45  cm.  A  Si 
photodetector  (Hammamatsu  S4751,  750  MHz  bandwidth)  measures  the  intensity 
of  light  emerging  from  the  interferometer.  The  sensitive  area  of  the  photodiode 
is  much  smaller  than  the  cross  sectional  area  of  the  laser  beam  so  only  a  fraction 
of  the  interferometer’s  output  is  detected.  The  size  of  the  detector  is  actually  an 
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Figure  3.2:  Photograph  of  the  experimental  apparatus.  The  black  box  left  of 
center  is  the  laser  housing.  The  path  of  the  laser  beam  is  highlighted  by  the  thick 
white  line.  The  coil  of  coax  cable  to  the  right  of  center  is  the  delay  line. 
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F  igure  3.3:  Schematic  diagram  of  the  electronic  feedback  loop  including  the  photo¬ 
diode  and  laser.  The  components  labelled  A-I  are  as  follows:  A  -  Hamamatsu  silicon 
photodiode  S4751;  B  -  MiniCircuits  directional  coupler  ZFDC-20-4;  C  -  MiniCir- 
cuits  amplifier  ZFL-1000L;  D  -  MiniCircuits  amplifier  ZFL-1000GH;  E  -  Coaxial 
cable  RG  58/U;  F  -  MiniCircuits  power  combiner  ZFSC-2-1W-75;  G  -  bias-T  in 
Thorlabs  laser  mount  TCLDM9;  H  -  Thorlabs  laser  diode  controller  LDC-500;  I  - 
Hitachi  laser  diode  HL6501MG. 

advantage  since,  due  to  aberrations  and  finite  precision  in  alignment,  more  than  one 
fringe  appears  within  the  beam  cross  section.  The  limited  aperture  allows  a  more 
spatially  uniform  held  to  fall  on  the  detector,  producing  better  fringe  visibility.  A 
neutral  density  filter  is  used  to  reduce  the  optical  power  reaching  the  photodiodes 
to  avoid  saturation. 

The  electronic  feedback  loop  which  begins  with  the  photodiode  and  ends  back  at 
the  laser  is  shown  schematically  in  Fig.  3.3.  The  part  numbers  of  the  components 
are  given  in  the  figure  caption.  The  photodiode  (labelled  A  in  Fig.  3.3)  produces  a 
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current  proportional  to  the  optical  power  falling  on  its  surface,  which  is  converted 
into  a  current  using  a  50  D  resistor.  The  voltage  across  this  resistor  is  transmitted 
down  a  coaxial  cable  (RU  58).  One  percent  of  the  signal  power  is  split  off  by 
a  directional  coupler  (B  in  Fig.  3.3)  and  sent  to  a  fast  oscilloscope  or  spectrum 
analyzer  to  monitor  the  state  of  the  system.  The  signal  propagating  down  the  main 
line  next  reaches  a  low-noise,  fixed-gain  amplifier  (C  in  Fig.  3.3),  a  DC-blocking 
capacitor  (220  pF),  a  variable  gain  amplifier  (D  in  Fig.  3.3),  and  a  second  DC- 
blocking  capacitor  (470  pF).  The  capacitors  serve  as  high  pass  filters  to  reduce  the 
loop  gain  at  frequencies  below  10  MHz  where  a  thermal  effect  enhances  the  laser’s 
sensitivity  to  frequency  modulation  .  A  power  splitter/combiner  (labelled  F  in  3.3) 
combines  the  feedback  signal  with  any  external  modulation  or  control  signal  that  is 
present  and  the  resulting  voltage  is  applied  to  the  bias-T  input  on  the  laser  mount. 
The  bias-T  (G  in  Fig.  3.3)  converts  the  signal  into  a  current  and  adds  it  to  a 
DC  component  from  a  commercial  laser  driver  (H  in  Fig.  3.3)  to  form  the  laser’s 
injection  current. 

The  entire  system  is  fixed  on  an  optical  table  using  short  (2  inch)  mounts  for 
mechanical  stability.  This  stability  is  extremely  important  as  variation  in  the  path 
length  on  the  order  of  the  the  wavelength  of  the  laser  light  (0.65  /mi)  produces 
significant  power  variations  at  the  output  of  the  interferometer.  Furthermore,  the 
entire  apparatus  is  covered  by  an  insulating  box  to  reduce  thermal  expansion  or 
contraction  of  the  mirror  mounts  due  to  air  currents. 

3.2  Mathematical  Model 

In  this  section,  I  develop  a  mathematical  model  of  the  active  laser  interferometer 
with  delayed  bandpass  feedback.  The  device  consists  of  three  main  components:  a 
semiconductor  laser  that  provides  linear  electro-optical  conversion,  a  Mach-Zehnder 
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interferometer  that  provides  nonlinearity,  and  an  electronic  feedback  loop  that  pro¬ 
vides  time  delay  and  bandpass  filtering.  I  consider  each  of  these  components  in 
a  separate  subsection.  A  block  diagram  appears  at  the  top  of  each  subsection  to 
highlight  the  component  to  be  discussed.  I  combine  the  component  models  to  arrive 
at  a  mathematical  model  of  the  full  experimental  system.  I  do  not  include  external 
modulation  in  this  section.  In  Sec.  3.5,  the  model  is  extended  to  include  modulation 
by  a  driving  voltage. 

3.2.1  The  Semiconductor  Laser 


F  igure  3.4:  The  first  component  of  the  active  interferometer  with  delayed  bandpass 
feedback:  The  semiconductor  laser.  The  laser  provides  linear  conversion  of  changes 
in  the  injection  current  to  changes  in  the  optical  power  and  frequency. 

The  light  source  in  the  active  laser  interferometer  is  a  semiconductor  laser.  In¬ 
vented  in  1962  [97]-[100],  the  semiconductor  laser  (or  laser  diode)  has  since  become 
a  key  component  of  optical  telecommunications  systems.  Reasons  for  its  suitability 
include  its  compact  size  (relative  to  other  lasers)  and  the  fact  that  its  power  can  be 
directly  modulated  through  variation  of  the  injection  current  [95].  These  features 
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also  make  it  suitable  for  my  purposes. 

Most  types  of  laser  involve  three  basic  elements:  a  pumping  mechanism,  a 
gain  medium,  and  a  resonator.  The  pumping  mechanism  adds  energy  to  the  gain 
medium,  maintaining  it  in  an  excited  state  (he.,  a  population  inversion).  Photons 
propagating  through  the  excited  medium  stimulate  coherent  emission  of  more  pho¬ 
tons.  The  resonator  ensures  that  each  photon  passes  several  times  through  the  gain 
medium,  and  is  thereby  amplified  many  times  before  escaping.  Light  leaves  the 
resonator  in  the  form  of  an  intense,  coherent  beam. 

All  three  of  these  elements  are  present  in  a  diode  laser.  The  laser  is  essentially  a 
semiconductor  p-n  junction  (or  homo junction) .  Light  is  emitted  when  electrons  and 
holes  recombine  in  the  depletion  region  of  the  junction  which  constitutes  the  gain 
medium.  The  pumping  mechanism  is  an  injected  current  that  produces  electron- 
hole  pairs.  Two  cleaved  facets  of  the  semiconductor  crystal  form  a  Fabry-Perot 
resonator.  A  schematic  drawing  of  a  homojunction  diode  laser  is  shown  in  Fig. 
3.5(a).  Modern  laser  diodes  differ  from  simple  p-n  junctions  because  additional  lay¬ 
ers  of  n-type  or  p-type  are  added  (forming  a  heterojunction )  to  confine  the  injected 
current  to  the  active  region.  The  various  layers  each  have  a  different  index  of  refrac¬ 
tion  and  conveniently  constitute  a  dielectric  waveguide  that  confines  the  optical  held 
to  the  gain  region.  The  entire  sandwich  of  semiconductor  layers  is  referred  to  as  a 
heterostructure.  One  type  of  GaAlAs  heterostructure  is  shown  in  Fig.  3.5(b).  If  the 
width  of  the  active  layer  is  comparable  the  DeBroglie  wavelength  of  the  electrons 
(~  10  nm),  carrier  motion  normal  to  the  active  layer  is  confined  to  a  quantum  well. 
So  called  quantum  well  lasers  tend  to  have  high  gain  at  lower  current  densities.  The 
diode  used  in  my  interferometer  is  a  multi-quantum  well  laser. 

The  dynamical  behavior  of  a  semiconductor  laser  can  be  modelled  by  rate  equa¬ 
tions,  a  set  of  first-order  differential  equations  that  describe  the  evolution  of  the 
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Figure  3.5:  (a)  Simplified  schematic  of  a  semiconductor  laser,  (b)  Schematic 
of  a  typical  “buried  heterostructure”  laser.  Cladding  layers  provide  both  carrier 
confinement  and  lateral  optical  waveguiding. 


34 


optical  field  and  carrier  density.  The  dynamics  of  a  single-mode  double  heterostruc¬ 
ture  semiconductor  laser  is  described  by  [95,  96] 

S  =  {G-'r)S  +  Rap,  (3.2) 

N  =  -  —  7eiV  —  GS,  (3.3) 

Q 

4>  =  —  —  (cu0  —  hi)  +  ~Pc  (G  —  7) ,  (3.4) 

I1  a  z 

where  S  and  N  are  the  number  of  photons  and  carriers,  respectively,  in  the  active 
layer,  (f)  is  the  phase  of  the  optical  held,  and  I  is  the  injection  current.  The  gain  G 
is  a  function  of  both  the  photon  and  carrier  numbers,  and  the  spontaneous  emission 
rate  Rsp  is  a  function  of  the  carrier  number.  The  the  loss  due  to  absorption  or  trans¬ 
mission  through  the  reflective  facets  is  denoted  by  7,  q  is  the  charge  of  an  electron, 
7e  is  the  rate  of  non-radiative  recombination  of  electron-hole  pairs,  jl  is  the  effective 
index  of  refraction  for  the  lasing  Fabry-Perot  mode,  ng  is  the  group  refractive  index, 
ojq/2ti  is  the  steady  state  optical  frequency,  12  is  the  cavity  resonance  frequency,  and 
(3C  is  the  linewidth  enhancement  factor.  The  optical  frequency  of  the  laser  is 

v  =  -^{uo  +  (j>)-  (3-5) 

Equations  3. 2-3. 4  constitute  a  detailed  model  of  the  semiconductor  laser.  How¬ 
ever,  they  contain  numerous  parameters  that  are  difficult  to  determine  empirically 
from  a  comercially  packaged  laser.  Also,  they  do  not  accurately  describe  the  effects 
of  injection  current  modulation  because  they  do  not  include  package  and  chip-related 
impedances  that  tend  to  shunt  current  around  the  active  layer  [94],  Fortunately, 
the  full  rate  equation  model  is  not  always  necessary.  When  the  injection  current  I 
is  fixed,  the  rate  equations  display  a  steady  state.  Below  a  threshold  current  value 
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Ith,  the  steady  state  power  is  nearly  zero  (except  for  a  small  contribution  from  spon¬ 
taneous  emission).  Above  threshold,  the  steady  state  power  is  linearly  proportional 
to  the  injection  current.  The  laser  approaches  the  steady  state  from  an  arbitrary 
initial  condition  via  damped  oscillations  with  a  damping  rate  Tr.  The  frequency  of 
these  oscillations,  denoted  by  Or,  is  known  as  the  relaxation  oscillation  frequency. 
Typically,  Or/27t  ~  Tr  ~  2  GHz.  A  simpler  model  of  the  semiconductor  laser  is 
applicable  when  the  injection  current  is  modulated  at  frequencies  well  below  Or,  as 
described  below. 

Agrawal  and  Dutta  [95]  show  that  the  optical  power  and  phase  respond  sinu¬ 
soidally  when  the  injection  current  is  weakly  modulated  around  a  constant  value 
such  that  I(t)  =  I dc  +  SI  sin (umt).  The  amplitude  6S  and  phase  6$  of  the  photon 
number  fluctuation  are  given  by 


6q  =  tan  1 


2T  Ru>m 

Ws  +  r|) 


\  \  (^R  +rR  )  J  ) 

where  S  is  the  steady  state  power  and  Gn  =  d G/dN.  Note  that  these  expressions 


reduce  to 
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Ss~  +  V  100  MHz  y  2ir' 

in  the  limit  where  ojm  «  \JG?r  +  Thus,  if  the  modulation  frequency  is 
small  enough  compared  to  the  relaxation  oscillation  frequency,  the  laser’s  opti¬ 


100  MHz  J  2v r 


cal  power  varies  linearly  with  the  injection  current  and  with  no  phase  lag.  A 
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similar  result  holds  for  the  optical  phase.  Agrawal  and  Dutta  give  the  phase  re¬ 
sponse  to  current  modulation  in  terms  of  the  resulting  change  in  optical  frequency 
6v{t)  =  5^0  sin  (umt  +  9C),  where 


6un  =  61 


PcG 


N 


Anq 


\ 


ut 


+  r2P 


Wh  +  r  If 


1  - 


+rR 
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4FR 


Op  +Fr  )  (s^r  +rj( ) 


(3.10) 


6C  =  tan  1  (  ^  )  +  tan  1 


where 
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2r  Rujm 


(fi|  +  r|) 
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Tp  =  ^r~sjp~2  GHz- 


(3,11) 


(3.12) 


In  the  limit  where  um  «  \/G2R  +  these  equations  reduce  approximately  to 


8v n  = 


PcG. 
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dvr q  (^|  +  r|) 
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(3.13) 
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rad 
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(3.14) 


rr.  (n|  +  r DJ~'"  v^’noo MHzy  2x 

Thus,  if  the  current  is  modulated  slowly  enough,  the  laser  may  be  viewed  as  a  simple 
linear  current-power  and  current-frequency  transducer. 

To  ensure  that  the  system  remains  in  the  regime  of  slow  modulation,  I  determine 
the  relaxation  oscillation  frequency  experimentally.  I  take  advantage  of  the  well 
known  fact  that  the  power  spectrum  of  the  intensity  noise  of  a  diode  laser  is  strongly 
peaked  around  the  relaxation  oscillation  frequency  [95].  With  the  laser  emitting  in 
a  steady  state,  the  intensity  noise  is  detected  by  a  fast  photodiode  (New  Focus  1580- 
B,  12  GHz  bandwidth).  The  frequency  content  of  the  signal  from  the  detector  is 
visualized  using  a  high-speed  spectrum  analyzer  (Hewlett-Packard  8566B,  22  GHz 
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Injection  Current  (mA) 

Figure  3.6:  Square  of  relaxation  oscillation  frequency  as  a  function  of  injection 
current.  By  operating  with  a  large  DC  injection  current,  I  ensure  that  the  laser  is 
never  modulated  at  frequencies  near  the  relaxation  oscillation. 

bandwidth).  The  location  of  the  peak  in  the  spectrum  is  recorded  as  the  injection 
current  of  the  laser  is  varied.  Figure  3.6  shows  the  result.  The  experimental  data 
reflect  the  well-established  relation  [95] 

oc  i-ith.  (3.15) 

Based  on  the  slope  obtained  from  a  least  squares  fit,  I  determine  the  relaxation 
oscillation  frequency  at  the  operating  current  of  74  mA  used  in  the  experienrents 
is  2.7  GHz,  much  higher  than  the  bandwidth  of  several  of  the  components  in  the 
electronic  feedback  loop.  Thus,  it  may  be  safely  assumed  that  the  laser  is  never 
modulated  at  a  rate  close  to  the  relaxation  oscillation  frequency. 

Based  on  the  results  of  the  preceeding  paragraphs,  I  approximate  the  laser’s 
behavior  as 

P  =  cr  (I  —  Ith)  =  a  (IAC  +  I dc  —  hh) ,  (3.16) 

v  —  V  (I  ~  hh)  =  V  (I ac  +  I dc  ~  hh) ,  (3-17) 
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where  P  is  the  total  optical  power  emitted  by  the  laser  diode,  and  cr  and  rj  are 
empirically  determined  constants.  The  injection  current  I  is  the  sum  of  a  constant 
bias  current  I^c  and  a  small  fluctuating  current  I  ac-  Fc>r  convenience,  the  constant 
terms  can  be  lumped  to  obtain 

P  =  vIac  +  Pqi  (3.18) 

V  =  7]  I  ac  +  v0,  (3.19) 

where  P0  =  a  (I DC  ~  hh)  and  v o  =  V  {I DC  ~  hh)-  Since  both  quantities  are  linear 
functions  of  I^c,  the  frequency  can  be  expressed  as  a  function  of  the  output  power 

zy  =  1  (p  -  p0)  +  u0.  (3.20) 

cr 

Equations  3.18  and  3.20  constitute  a  more  practical  model  of  the  semiconductor 
laser  because  they  contain  only  a  few  parameters,  all  of  which  are  experimentally 
accesible.  This  model  inherits  two  weakness  from  the  rate  equation  model  from 
which  it  was  derived.  First,  the  rate  equation  model  only  includes  one  mechanism 
for  frequency  modulation,  i.e.  variation  of  the  carrier  number.  In  reality,  there  is 
a  second  physical  process  that  produces  frequency  modulation  when  the  injection 
current  is  varied.  That  process  is  a  thermal  effect  that  enhances  frequency  mod¬ 
ulation  below  10  MHz  by  as  much  as  a  few  orders  of  magnitude  [91,  92,  96].  As 

mentioned  earlier,  the  capacitors  in  the  feedback  loop  in  the  active  laser  interferom¬ 
eter  introduces  greater  loss  at  low  frequencies  to  counteract  this  effect.  Therefore, 
I  do  not  account  for  the  thermal  effect  explicitly  in  my  model.  However,  care  must 
be  taken  in  measuring  the  constant  rj  as  it  will  have  a  different  value  below  10  MHz 
than  above. 

The  second  weakness  of  this  model  is  that  it  neglects  the  important  physical 
phenomenon  of  phase  noise.  Due  to  spontaneous  recombination  in  the  active  region, 
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the  phase  of  the  electric  held  is  subject  to  discontinuous  jumps  [95].  As  a  result, 
the  frequency  of  the  laser  fluctuates  somewhat  around  the  value  predicted  by  the 
deterministic  model.  The  effect  of  the  interferometer  on  this  phase  noise  is  discussed 
in  the  next  subsection. 

3.2.2  The  M  ach-Zehnder  I  interferometer 


F  igure  3.7:  The  second  component  of  the  active  interferometer  with  delayed  band¬ 
pass  feedback:  The  Mach-Zehnder  interferometer.  The  interferometer  provides  the 
nonlinearity  needed  to  obtain  chaos. 

The  coherent  light  waves  light  emitted  by  the  diode  laser  are  directed  into  an 
interferometer.  As  will  be  shown  below,  the  interferometer  outputs  a  beam  whose 
intensity  is  a  nonlinear  function  of  the  input  optical  power  and  frequency.  Thus,  this 
component  provides  nonlinearity  to  the  system,  a  necessary  property  for  producing 
complex  dynamics.  In  this  subsection,  I  explain  the  principle  of  operation  of  the 
Mach-Zehnder  interferometer  in  detail. 

A  Mach-Zehnder  interferometer  takes  a  single  laser  beam  as  input,  splits  it  into 
equal  intensity  beams  that  travel  in  separate  arms,  and  recombines  them  after  they 
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travel  different  distances  (see  Fig.  3.8(a)).  To  model  this  process,  I  assume  the 
input  held  is  of  the  form 

Ein(x,  y,  z,  t )  =  E(y,  z)e~i^-kx\  (3.21) 

with  optical  frequency  v  =  o;/27t  and  intensity 

ha  =  ^E2,  (3.22) 

where  c  is  the  speed  of  light  and  eo  is  the  permittivity  of  free  space.  Since  there  are 
two  50-50  beam  splitters  (labelled  BS1  and  BS2  in  Fig.  3.8(a))  in  the  system,  there 
are  actually  four  separate  beam  paths,  as  shown  in  Fig.  3.8(b).  Two  of  these  paths 
meet  at  each  output.  Reflection  from  each  air-to-glass  interface  adds  a  tt  phase  shift 
to  the  electric  held.  Note  that  the  reflection  in  path  4  at  BS2  is  at  a  glass-to-air 
interface  so  no  phase  shift  occurs.  Paths  1  and  3  are  longer  than  paths  2  and  4  by  a 
distance  l  =  cA  where  A  is  the  difference  in  propagation  time  through  the  two  arms 
of  the  interferometer.  The  electric  helds  following  each  path  can  then  be  expressed 
as 

Ei  (x, y,  z,  t)  =  — 1  E(y,  (3.23) 

E,(x,y,z,t)  =  -i  E(y,z)e-i^~k‘\  (3.24) 

E3(x,  y,  z,  t)  =  i  E(y,  *)e-fM-‘<*-‘»,  (3.25) 

Ea(x,  y,  z,  t)  =  ~E(y,  (3.26) 

where,  Et  is  the  held  that  follows  path  i.  The  intensity  at  either  output  is  the 
coherent  combination  of  the  helds  associated  with  the  two  paths  that  reach  it. 
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Fields  Ei  and  E2  are  combined  at  output  1  so  the  intensity  at  output  1  is 


i1  =  c^\e1  +  e2\2, 


_  ce0 
2 


=  -^-E2  [1  +  cos (kl)] , 
=  Y  i1  +  cos(cnA)] , 


1  1 


_  i  _  /  —i(cot—k(x—l))  i(ujt—kx)  .  —i(u>t—kx)i(ujt—k(x—l))\ 

2  +  4  ve  e  “*■ e  e  7  ’ 


(3.27) 

(3.28) 

(3.29) 

(3.30) 


where  I  have  used  k  =  cu  and  kl  =  k(cA )  =  uA.  The  fields  E3  and  E^  are  combined 
at  output  2  giving  an  intensity 


h  —  — ~  \Es  +  E4\2 


=  ceo  ^2 
2 


=  [!  -  cos  (kl)] , 

=  Y  t1  -  cos(cnA)] . 


1 

4 


_  _  _  ^—i(ujt—k(x—l))^i(Lut—kx) 


|  ^—i(ujt—kx)^i(ujt—k(<x—l))'j 


(3.31) 

(3.32) 

(3.33) 

(3.34) 


This  model  of  the  interferometer  can  be  connected  to  the  model  of  the  diode 
laser  since  it  is  the  source  of  the  input  intensity  Iin.  In  the  laser  model,  the  output 
light  is  described  in  terms  of  the  optical  power  emitted.  The  optical  power  P  and 
intensity  /  ( x ,  y )  of  the  light  emitted  by  the  laser  are  related  by 


P  = 


(3.35) 


where  the  integration  is  performed  over  the  elliptical  cross  section  of  the  laser  beam. 
The  important  feature  of  Eq.  3.35  is  the  proportionality  between  P  and  /.  This 
result,  together  with  Eqs.  3.30  and  3.34,  implies 

Pi  =  P  [1  +  cos(n;A)] ,  (3.36) 
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(3.37) 


P2  =  yP  [1  -  COs(cuA)]  , 

where  I  have  included  the  factor  Q?  to  account  for  absorption  by  the  neutral  density 
hlter  between  the  laser  and  the  interferometer. 

Equation  3.20  relates  the  optical  frequency  v  of  the  laser  to  the  optical  power 
P.  Thus,  with  a  laser  diode  as  the  light  source,  the  outputs  of  the  interferometer 
are 

Cf 


Pi  =  yP  {1  +  cos  [<x(p-Po)  +  </>]}, 
P2  =  yP  {1  -  cos  [<x(P-Po)  +  </>]}, 


(3.38) 


where  a  =  2t\Ati/o  and  (j)  =  2'KAr)/av0.  The  constant  offset  <j)  of  the  argument  to 
the  cosine  can  be  tuned  by  varying  the  DC  component  of  the  injection  current.  I 
choose  to  set  this  value  to  m  ±  tt/2  where  m  is  an  even  integer.  If  (f)  =  — vr/2, 


Pi  =  ^P{l  +  sin[a(P-P0)]}, 
P2  =  yP{l-sin  [a  (P-P0)]}. 


(3.39) 

(3.40) 


Choosing  the  opposite  sign  of  (f)  simply  flips  the  sign  of  the  cosine  terms.  The 
constant  parameter  a  determines  the  sensitivity  of  the  interferometer  and  can  be 
tuned  by  varying  the  path  imbalance  A.  With  a  large  value  of  a,  the  interferometer 
will  amplify  small  fluctuations. 

Up  to  this  point,  I  have  assumed  the  light  waves  in  the  interferometer  are  per¬ 
fect  plane  waves  and  that  the  beamsplitters  divide  them  exactly  in  half.  In  reality, 
neither  of  these  conditions  are  strictly  true.  Curvature  of  the  wavefronts  produces 
variation  in  the  degree  of  interference  across  the  beamwidth.  Typically,  the  inter¬ 
ference  pattern  is  a  set  of  curved  stripes.  Careful  alignment  of  the  interferometer 
makes  the  stripes  thicken  until  only  a  couple  are  visible  over  the  beam  width.  How¬ 
ever,  the  ideal  of  a  uniform  held  is  not  experimentally  attainable.  Likewise,  since 
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the  beamsplitters  are  not  perfect  50-50  splitters,  the  waves  in  each  arm  have  differ¬ 
ent  amplitudes  so  one  wave  cannot  completely  cancel  out  each  the  other.  Both  of 
these  effects  tend  to  reduce  the  visibility  b  of  the  fringes  defined  as 


b  = 


p  . 

±  mir 


+  Pu 


(3.41) 


I  account  for  this  effect  by  adding  a  parameter  to  Eqs.  3.39  and  3.40  that  reduces 
the  visibility.  Specifically, 


^-P{l  +  bsm[a  (P-P0)]}, 

(3.42) 

^P{l-6sin  [a(P-Po)]}. 

(3.43) 

Equations  3.42,  3.43,  and  3.18  consitute  a  practical  model  of  the  laser  and  in¬ 
terferometer.  Before  proceeding  to  model  the  feedback  loop,  several  caveats  must 
be  noted.  First,  since  the  interferometer  sensitivity  a  is  directly  proportional  to  the 
path  imbalance,  at  first  glance  it  appears  that  a  can  be  set  to  an  arbitrarily  large 
value.  Unfortunately,  the  interferometer  also  amplifies  phase  noise  in  the  laser  field. 
To  see  why,  assume  that  the  phase  noise  can  be  modelled  by  adding  a  stochastic 
term  £(i)  to  Eq.  3.20  giving 

v=1±[(P-P0)  +  v0  +  at)}.  (3.44) 

a 

The  output  of  the  interferometer  is  then 

P1  =  ^-P{  1  +  cos  [a  ( P  -  P0)  +  a£(t)  +  </>]}.  (3.45) 

Increasing  the  sensitivity  of  the  interferometer  proportionally  increases  the  effect  of 
phase  noise.  Thus  there  is  a  tradeoff  between  sensitivity  and  noise.  Experimentally, 
I  chose  to  set  the  sensitivity  to  a  value  such  that  amplitude  of  noise  fluctuations  is 
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roughly  less  than  10%  of  Pq.  This  amount  of  noise  may  be  large  enough  to  have  a 
significant  effect  on  the  dynamics  of  the  system. 

Second,  I  have  assumed  that  the  light  waves  reaching  one  output  through  two 
different  paths  have  the  same  optical  frequency.  This  assumption  is  reliable  as  long 
as  the  optical  frequency  is  not  modulated  significantly  in  a  period  of  time  compara¬ 
ble  to  A,  the  difference  in  propagation  times  through  the  two  paths  being  combined. 
The  path  difference  in  my  interferometer  is  ~  45  cm  making  the  difference  in  prop¬ 
agation  times  ~  1.5  ns.  Oscillations  in  the  system  at  frequencies  comparable  to  or 
greater  than  1/(1. 5  ns)  or  670  MHz  cannot  be  explained  by  this  model.  In  order  to 
scale  this  device  to  higher  frequencies,  either  the  path  difference  must  be  shortened 
or  a  new  model  must  be  developed. 

Third,  due  to  aberrations  in  the  wavefronts  and  finite  precision  in  aligning  the 
optical  components  making  up  the  interferometer,  the  width  of  the  interference 
fringes  is  less  then  the  width  of  the  laser  beam  itself,  resulting  in  some  spatial  varia¬ 
tion  in  intensity  over  a  cross  section  of  the  beam  emerging  from  the  interferometer. 
Typically,  about  two  fringes  fit  within  the  beam  width.  However,  the  output  of  the 
interferometer  is  detected  by  a  photodiode  whose  active  region  of  is  much  smaller 
than  the  beam  width  so  only  a  single  fringe  is  observed.  As  a  result,  since  much  of 
the  beam  is  clipped,  only  a  fraction  of  the  power  is  detected. 

Fourth,  many  types  of  interferometer  exist,  e.g.  the  well-known  Michealson 
interferometer,  that  could  be  used  as  the  nonlinear  element  in  this  system.  I  chose 
to  use  a  Mach-Zehnder  interferometer  because  the  components  are  arranged  in  such 
a  way  that  no  light  is  reflected  back  into  the  laser  cavity.  Such  back  reflections  are 
a  common  cause  of  instability  in  diode  lasers  [80].  Ohtsubo  et  al.  [43]  investigated 
a  similar  interferometer-based  device  but  used  a  Twyman-Green  interferometer  to 
provide  nonlinearity.  As  a  result,  they  had  to  include  an  expensive  optical  isolator 
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in  order  to  guard  against  instabilities  induced  by  back  reflections. 

The  light  exiting  output  1  of  the  interferometer  falls  on  a  photodiode  that  con¬ 
verts  the  optical  signal  into  an  electronic  signal.  This  signal  propagates  through  the 
electronic  feedback  loop  which  is  the  subject  of  the  next  subsection. 

3.2.3  The  Electronic  Feedback  Loop 


Figure  3.9:  The  third  component  of  the  active  interferometer  with  delayed  band¬ 
pass  feedback:  The  electronic  feedback  loop.  Time  delay  and  bandpass  filtering  are 
produced  by  the  feedback  loop. 

The  electronic  feedback  loop  converts  the  optical  intensity  emitted  by  the  in¬ 
terferometer  into  an  electrical  signal  and,  after  a  delay  period,  modulates  that  sig¬ 
nal  onto  the  injection  current  of  the  diode  laser.  The  loop  contains  several  linear 
elements  with  amplification,  attenuation,  and  filtering  characteristics.  High-pass 
filtering  is  provided  by  the  two  capacitors  and  the  bias-T.  Low-pass  filtering  is  pro¬ 
vided  by  the  limited  bandwidth  of  the  photodiode,  diode  laser,  and  the  electrical 
connections  to  the  laser.  As  was  the  case  with  semiconductor  laser,  developing 
a  first  principles  model  of  the  feedback  loop  requires  significant  complexity  and 
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Figure  3.10:  Block  diagram  of  feedback  loop.  Light  from  the  interferometer  is 
converted  into  an  electrical  signal  by  the  photodiode  (PD).  All  electrical  components 
are  lumped  into  one  amplifier  with  gain  7,  one  low-pass  filter  with  time  constant 
Tl,  one  high-pass  filter  with  time  constant  th,  and  a  time  delay  of  td-  Finally,  the 
electrical  signal  is  converted  back  into  light  by  the  laser  diode  (LD). 

knowledge  of  parameters  that  are  difficult  to  obtain  empirically  from  commercial 
devices  such  as  the  bias-T  and  laser  mount.  Therefore,  I  will  again  develop  a  highly 
simplified  but  practical  model.  I  will  assume  the  amplification,  attenuation,  and 
filtering  of  the  various  of  the  components  of  the  feedback  loop  can  be  lumped  into 
one  single-pole  low-pass  Liter,  one-single  pole  high-pass  Liter,  and  one  ampliLer. 

Figure  3.10  shows  a  block  diagram  of  the  lumped  model.  Light  from  the  inter¬ 
ferometer  falls  on  the  detector  (PD  in  Fig.  3.10)  producing  a  current  proportional 
to  the  incident  intensity.  Based  on  the  interferometer  response  given  in  Eq.  3.42, 
the  photocurrent  is 

Met (t)  =  —P(t)  {1  +  6 sin  [a  ( P(t )  -  PQ)]}  ,  (3.46) 

where  p  is  the  sensitivity  of  the  photodiode  and  £ c  is  the  fraction  of  the  power  in  the 
laser  beam  that  actually  falls  on  the  detector.  The  photocurrent  is  converted  into 
a  voltage  by  a  resistance  r  =50  Q  and  amplihed  by  a  factor  £  producing  a  voltage 

Vdet  =  {1  +  frsin  [a  (P  -  P0)]}  ,  (3.47) 


48 


where 


pCcCf^t  1 
2  y/2' 


(3.48) 


The  extra  factor  of  2~l/-  is  included  to  account  for  losses  in  the  directional  coupler 
(labelled  B  in  3.3)  and  power  splitter /combiner  (labelled  F  in  3.3),  respectively. 
This  voltage  delayed  for  a  length  of  time  Td  before  reaching  the  low-pass  filter.  In 
general,  a  single  pole  low-pass  filter  is  described  by  the  equation 


tlV  =  -V  +  Vin(t), 


(3.49) 


where  Vin  is  the  input  to  the  filter,  V(t)  is  the  output,  and  ry  is  the  response  time. 
In  this  case,  the  input  is  Vc\et(t  —  7y>)  so  the  output  of  the  low-pass  filter  is  given  by 
the  equation 

tlV  =  -V  4-  jPm  {1  +  b  sin  [a  (P^  -  P0)]}  ,  (3-50) 

where  PTD  =  Pit—To)-  The  output  of  the  low-pass  hlter  is  the  input  to  the  high-pass 
filter.  A  single  pole  high-pass  hlter  is  described  by  the  equation 

U  = - b  Uin,  (3.51) 

th 

where  Uin  is  the  input  to  the  hlter,  U (t)  is  the  output  voltage,  and  Th  is  the  response 
time.  In  this  case,  Uin  =  V  so 

U  =  — —  +  V.  (3.52) 

th 

I  assume  this  voltage  is  proportional  to  the  current  Iac  appearing  in  Eqs.  3.18  and 
3.19  so  that 


P(t)  =  aRU(t )  +  P0, 


(3.53) 
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where  R  is  the  constant  of  proportionality  and  has  dimensions  of  impedance.  Taking 
the  derivative  with  respect  to  time  of  this  equation  and  substituting  Eq.  3.52  for  U 
I  obtain 

P=-—{P-  P0)  +  kV,  (3.54) 

Tff 

where  k  =  aR.  Equation  3.54  describes  the  evolution  of  the  optical  power  output 
by  the  laser  diode.  The  light  from  the  laser  travels  through  the  interferometer  and 
ultimately  reaches  the  detector,  thereby  closing  the  feedback  loop. 

3.3  Experimental  Determination  of  M  odel  Param¬ 
eters 


Equations  3.50  and  3.54  together  constitute  the  complete  model  of  the  active  in¬ 
terferometer  with  time-delayed  bandpass-filtered  feedback.  For  reference,  I  restate 
them  here 

P  =  -—{P-Pq)  +  kV, 
th 

tlV  =  -V+jPTD{l+bsm[a(PTD-P0)}}.  (3.55) 


The  model  contains  ten  parameters.  Table  3.1  lists  their  values  as  determined 
experimentally.  In  this  section,  I  describe  the  measurements  I  perform  to  obtain 
these  values. 

The  feedback  loop  gain  7  depends  on  a  number  of  parameters  that  are  straight 
forward  to  measure.  Specifically,  according  to  Eq.  3.48, 


(50  ohms)  p(c(f£, 
2V2 


(3.56) 


where  p  is  the  senitivity  of  the  photodiode,  Qc  is  the  fraction  of  the  laser  beam  that 
falls  on  the  small  sensitive  region  of  the  photodiode,  (p  is  the  fraction  of  light  not 
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Parameter 

Description 

Value 

Tl 

Low-pass  filter  time  constant 

0.66  ±0.05  ns 

Th 

High-pass  filter  time  constant 

22  ±0.5  ns 

Td 

Feedback  delay  time 

19.1  ±0.1ns 

K 

Voltage  to  power  conversion  at  bias-T 

4.8  ±0.1  x  10-3  mW/mV 

a 

Interferometer  sensitivity 

1.89  ±0.05mW_1 

Po 

Optical  power  operating  point 

26  ±0.5  rnW 

b 

Fringe  visibility 

0.8  ±0.02 

P 

Photodiode  sensitivity 

0.44  ±0.02  mA/mW 

Cc 

Fraction  of  total  optical  power  detected 

0.21  ±0.02 

C F 

Fraction  of  power  transmitted  by  NDF 

0.47  ±0.02 

Table  3.1:  Physically  relevant  values  of  model  parameters 


absorbed  by  the  nuetral  density  filter,  and  £  is  the  voltage  gain  of  the  amplifiers. 
The  amplifier  gain  £  is  determined  by  input  a  known  radio- frequecy  voltage  and 
observing  the  output.  The  gain  has  some  frequency  dependence  (less  than  3  dB) 
within  the  bandwidth  (0.01-1  GHz)  so  I  use  specifically  the  value  at  52  MHz,  the 
fundamental  frequency  of  oscillation  reported  in  the  next  section. 

The  photodiode  sensitivity  p,  and  losses  due  to  beam  clipping  and  the  neutral 
density  filter  are  determined  by  applying  a  known  amount  of  light  and  observing  the 
resulting  photocurrent  as  follows.  The  light  source  for  these  measurements  is  the 
laser  diode.  The  total  power  out  with  an  injection  current  of  74  rnA  is  determined 
to  be  Po  =  26  mA  using  a  large-area  photoreceiver  (New  Focus  2031).  I  place  the 
neutral  density  filter  between  the  laser  and  the  photoreceiver  to  find  the  fraction 
of  power  transmitted,  I  then  replace  the  photoreceiver  with  the  photodiode.  I 
add  a  lens  to  the  beam  path  to  obtain  a  smaller  beam  width  at  the  photodiode. 
The  focusing  ensures  that  all  the  optical  power  in  the  beam  reaches  the  sensitive 
area  of  the  photodiode.  However,  the  photodiode  is  not  positioned  precisely  at  the 
focal  point  of  the  lens  to  prevent  localized  saturation  of  the  detector.  The  current 
produced  by  the  photodiode  flows  through  a  50  D  resistance  and  the  resulting 
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Figure  3.11:  Setup  for  open  loop  measurements  used  to  determine  the  transfer 
characteristics  of  the  feedback  loop. 

voltage  is  recorded  to  determine  p.  Next,  I  place  the  laser  at  the  input  to  the 
interferometer  and  the  photodiode  at  the  output.  The  photocurrent  is  measured 
one  more  time  to  determine  £ c ■ 

The  parameters  tl,  th,  and  k  determine  the  pass  band  of  the  electronic  feedback 
loop.  I  obtain  reasonable  values  for  these  parameters  by  measuring  the  open  loop 
transfer  characteristic  of  the  feedback  loop.  The  loop  is  broken  by  disconnecting 
the  photodetector  and  repositioning  it  directly  in  front  of  the  laser  as  shown  in 
Fig.  3.11.  A  signal  generator  is  connected  in  place  of  the  photodetector  so  that 
a  sinusoidal  voltage  of  known  amplitude  and  frequency  is  injected  into  the  opened 
feedback  loop  and  eventually  to  the  laser.  The  output  of  the  laser  is  measured 
immediately  before  the  interferometer  using  the  photodetector  previously  removed 
from  the  loop.  This  open  loop  configuration  is  described  by  a  simple  modification 
of  Eq.  3.55.  In  this  case,  the  input  to  the  low-pass  filter  is  a  sinusoidal  voltage  of 
the  form  V%n  =  Voe*nt,  where  Vq  and  Q  are  the  amplitude  and  frequency  of  the  input 


52 


voltage,  respectively.  The  voltage  that  modulates  the  laser  is  then  given  by 


(3.57) 


rLV  =  -V  +  —=£V0eint. 

v  2 

Assuming  a  sinusoidal  response,  the  solution  to  this  equation  is 

V(t)  =  eiQt • 

1  +  iVLtl 

Inserting  this  result  into  Eq.  3.55  and  assuming  the  response  of  the  optical  power 
is  of  the  form  P  =  PDq  +  Pac where  Pdc  is  a  constant  offset  and  Pac  is  a 
complex  amplitude,  I  obtain 


(3.58) 


ittPACeiCit  =  (i Pdc  +  PACeiat  -  P0)  + 

Th  1  +  lilTL 


Equating  the  constant  terms  I  find 


(3.59) 


Pdc  =  Pq. 


(3.60) 


Equating  the  oscillating  terms,  cancelling  common  factors,  and  solving  for  P,\c  gives 


Pac  = 


Qth 


(1  +  ifirL) 


(3.61) 


Finally,  solving  for  the  magnitude  of  PAc  and  rearranging  I  obtain 

\Pac\ 


K 


ISM 


q2th 


(3.62) 


(1  +  Wrl) 


This  equation  can  be  compared  with  the  experimentally  determined  open  loop  re¬ 
sponse  in  order  to  obtain  physically  relevant  values  of  the  parameters  k,  tl,  and  th- 
Figure  3.12  shows  both  the  measured  (circles)  and  theoretical  (dotted  line)  response 
of  the  open  loop  system  to  sinusoidal  driving  at  two  different  values  of  the  amplifier 
gain  £.  The  theoretical  curve  is  produced  using  the  values  k  =  4.8  x  ICE3  mW/rnV, 
tl  =  0.66  ns,  and  th  =  22  ns  chosen  by  varying  k,  tl,  and  th  hr  order  to  achieve 
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Figure  3.12:  Frequency  response  of  open  loop  system  at  two  different  gain  values: 
(a)  7  =  14.6,  (b)  7  =  4.8. 

a  good  fit.  With  these  values,  the  theoretical  curve  fits  to  within  a  few  percent  of 
the  experimental  data  everywhere  except  in  the  region  between  80  and  140  MHz. 
The  dip  in  the  response  here  is  also  apparent  in  the  response  of  the  laser  to  direct 
modulation  through  the  bias-T.  Therefore,  it  is  apparently  caused  by  either  the 
bias-T  or  the  electronics  inside  the  laser  mount.  I  choose  to  ignore  this  discrepancy 
because,  as  the  next  two  sections  show,  the  simple  model  successfully  reproduces 
most  features  of  the  observed  dynamics. 

With  the  photodiode  placed  back  at  the  output  of  the  interferometer,  the  vis¬ 
ibility  b  is  determined  by  observing  minima  and  maxima  in  the  photodetector’s 
output  current  as  the  DC  injection  current  is  varied.  From  these  measurements, 
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Figure  3.13:  Setup  for  open  loop  measurements  used  to  determine  parameter  a. 

the  minimum  and  maximum  optical  power  is  calculated  for  several  fringes  and  Eq. 
3.41  is  evaluated  to  find  b. 

A  second  open  loop  measurement  is  made  to  determine  the  parameter  a,  in  this 
case  with  the  photodiode  placed  after  the  interferometer  as  shown  in  Fig.  3.13.  I 
assume  the  optical  power  entering  the  interferometer  is  P(t)  =  Pq  +  Pac PQt  where 
Pac  is  given  by  Eq.  3.61.  Inserting  this  into  Eq.  3.42,  I  obtain  the  power  output 
by  the  interferometer 

Pi  =  Y  (Po  +  PACeiQt )  {1  +  bsm  (aPACem)}  .  (3.63) 

Assuming  aPAc  «  1,  which  can  be  assured  by  keeping  the  drive  voltage  V0  small, 
I  approximate  Pi  as 

Pi  =  Y  (P°  +  PAceiQt)  { 1  +  baPACeint}  .  (3.64) 

I  collect  terms  according  to  oscillation  frequency  to  get 

Pi  =  Y  {Po  +  (1  +  Poba )  PACem  +  baPlcei2Qt}  .  (3.65) 
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The  amplitude  of  the  component  at  the  drive  frequency  Q  is 


—  (1  +  P0ba)  PAC.  (3.66) 

The  current  from  the  photodiode  is  measured  and  analyzed  using  a  spetrum  ana¬ 
lyzer.  Equation  3.66  is  then  used  to  calculate  a. 

The  remaining  parameter  is  the  time  delay  7y>  Direct  measurement  of  this 
parameter  is  not  easy  because  it  includes  contributions  from  several  electronic  com¬ 
ponents,  each  of  which  introduces  some  phase  lag  in  addition  to  simple  propagation 
delay.  Therefore,  I  resort  to  an  indirect  method  that  makes  use  of  the  model  itself.  I 
close  the  feedback  loop  and  allow  the  system  to  oscillate  periodically.  The  feedback 
gain  7  is  adjusted  so  that  the  system  is  just  beyond  the  Hopf  bifurcation.  In  the 
next  section,  I  present  a  linear  stability  analysis  of  this  bifurcation  that  relates  the 
observed  frequency  to  the  time  delay.  I  assume  the  correct  value  of  td  is  that  which 
gives  the  same  frequency  of  oscillation  at  the  Hopf  bifurcation  as  is  observed  in  the 
experimental  system. 

3.4  Dynamics  of  Closed-Loop  System  without  M  od- 
ulation 

To  characterize  the  behavior  of  the  active  interferometer  with  time-delayed  bandpass 
feedback,  I  examine  the  bifurcations  that  occur  as  I  vary  the  loop  gain  7  from  a 
low  to  a  high  value.  In  both  the  experiment  and  the  model,  I  observe  sequence  of 
bifurcations  beginning  with  a  steady  state  and  ending  with  broadband  chaos.  In  this 
section,  I  first  discuss  the  bifurcation  from  steady  state  to  a  periodic  oscillation.  I 
present  a  linear  stability  analysis  to  explain  the  fundamental  frequency  of  oscillation. 
Then  I  discuss  the  later  transition  from  periodic  oscillation  to  broadband  chaos. 

To  examine  the  bifurcations  of  a  dynamical  system  without  delay,  it  is  often 
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useful  to  visualize  trajectories  in  phase  space.  Much  information  can  be  gleaned 
from  the  geometry  and  topolgy  of  such  a  representation.  However,  the  phase  space 
of  a  delay  system  has  an  infinite  number  of  dimensions  and  is  difficult  to  visualize. 
Therefore,  I  simply  examine  time  series  and  power  spectra  in  order  to  explore  the 
dynamics  of  my  system.  The  two  variables  in  the  model,  P  and  V,  are  not  easily 
accesible.  The  power  emitted  by  the  laser  P  could  be  measured  by  putting  a  beam 
splitter  between  the  laser  and  the  interferometer.  However,  it  is  preferable  to  keep 
the  free  space  propagation  time  as  small  as  possible  to  minimize  the  response  time 
of  the  active  interferometer.  The  voltage  V  is  only  accessible  by  opening  the  laser 
mount  itself.  An  experimentally  convenient  quantity  is  the  output  from  the  unused 
port  of  the  interferometer.  This  quantity  is  compared  with  the  output  of  the  model 
using  Eq.  3.43.  All  measurements  described  in  the  following  subsections  are  from 
this  second  port  using  a  detector  identical  with  the  one  in  the  feedback  loop.  I  record 
the  voltage  signal  from  this  detector  in  the  time  domain  using  a  high-speed  digital 
oscilloscope,  and  in  the  frequency  domain  using  a  spectrum  analyzer.  Time  series 
from  the  model  are  obtained  by  numerical  integration  using  an  Adam-Bashforth- 
Moulton  predictor-corrector  algorithm  [101]. 

3.4.1  Hopf  Bifurcation 

With  the  gain  below  a  critical  value,  both  the  experiment  and  model  display  a 
steady  state.  Figure  3.14  shows  a  time  series  of  the  optical  power  detected  at  the 
second  output  of  the  interferometer  while  the  system  is  in  the  steady  state.  The 
fluctuations  are  attributed  to  noise  in  the  loop,  primarily  phase  noise  in  the  diode 
laser.  The  existence  of  the  steady  state  in  the  model  is  easily  demonstrated  by 
setting  the  derivatives  in  Eqs.  3.50  and  3.54  to  zero  and  solving  for  P  and  V .  The 
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Figure  3.14:  Time  series  of  optical  power  detected  at  second  port  of  interferometer 
with  feedback  gain  7  =  4.0  mV/mW.  The  system  is  in  a  steady  state  with  fluctations 
due  to  phase  noise. 

result  is 


P* 


(3.67) 


V*  =  7P0,  (3.68) 

where  P*  and  V*  denote  the  steady  state  values  of  P  and  V.  Assuming  the  value 
of  P0  as  given  in  Table  3.1,  the  power  out  of  the  second  port,  given  by  Eq.  3.43,  is 
6.11  rnW,  in  good  agreement  with  the  data  in  Fig.  3.14. 

When  the  gain  in  the  model  is  increased  through  the  critical  value  7*  =  5.34 
mV/mW,  the  steady  state  is  replaced  by  a  periodic  oscillation.  As  shown  by  the 
open  circles  in  Figure  3.15(a),  the  amplitude  of  this  oscillation  grows  smoothly  from 
zero  as  the  gain  is  increased.  In  the  range  just  beyond  the  bifurcation  point,  the 
amplitude  scales  as  y/j,  as  shown  in  Fig  3.15(b),  suggesting  a  supercritical  Hopf 
bifurcation  [48].  A  similar  transition  occurs  in  the  experiment  where  the  critical 
value  of  7  is  5.1  ±  0.5  mV/mW.  The  filled  circles  in  Fig  3.15(a)  show  how  the 
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Figure  3.15:  Experimental  and  simulated  data  showing  how  the  amplitude  of  the 
periodic  oscillation  increases  with  the  feedback  gain  7.  (a)  The  filled  circles  show  the 
experimentally  measured  amplitude.  The  large  error  bars  are  due  to  the  influence  of 
phase  noise  near  the  transition.  The  open  circles  show  data  from  the  model,  (b)  A 
closer  look  at  the  model  data  near  the  bifurcation  point.  The  line  is  a  least  squares 
fit  to  the  points  just  beyond  the  bifurcation.  The  linear  scaling  of  the  square  of  the 
amplitude  with  the  feedback  gain  indicates  a  supercritical  Hopf  bifurcation. 
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Figure  3.16:  Time  series  from  the  (a)  experiment  and  (b)  model  of  the  output 
at  the  second  interferometer  port  in  periodic  regime.  The  feedback  gain  is  7  =  6.8 
mV/mW. 

amplitude  increases  with  the  feedback  gain.  The  rather  large  uncertainty  is  due 
to  the  fact  that  the  phase  noise  obscures  the  transition  by  exciting  the  system  to 
oscillate  at  the  Hopf  frequency  even  when  it  is  slightly  unstable.  Figure  3.16  shows 
time  series  of  the  periodic  oscillation  from  the  experiment  and  model  with  7  =  6.8 
mV/mW. 

The  dominant  frequency  of  the  oscillation  is  51.5  MHz  ±  1  MHz.  This  frequency 
is  very  nearly  equal  to  l/iy>  By  varying  td  in  both  the  model  and  experiment,  I 
find  that  this  relation  holds  over  a  wide,  but  limited  range  of  Td-  Experimentally, 
td  is  varied  by  adding  or  subtracting  fixed  lengths  of  coaxial  cable  to  the  feedback 
loop.  The  first  experimental  value  of  td  =  19.1  ns  is  determined  according  to  the 
procedure  described  in  the  preceeding  section.  The  other  values  are  determined 
relative  to  this  value  by  adding  or  subtracting  propagation  time  appropriate  to 
the  change  of  cable  length.  Figure  3.17  shows  the  relation  between  the  observed 
frequency  and  Td  in  model  and  experimental  time  series.  In  general,  it  appears  the 
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frequency  is  approximately  ti/td  where  n  is  an  integer  that  increases  with  td-  The 
experimental  data  points  coincide  with  the  model  data  points  except  at  td  =  65.9 
ns.  Here  the  experimental  system  remains  on  the  n  =  2  branch  while  the  model  has 
switched  to  the  n  =  3  branch.  This  may  indicate  multistability  above  the  bifurcation 
point.  The  set  of  Hopf  frequencies  observed  represents  a  distinct  difference  between 
this  dynamical  system  and  those  like  the  Ikeda  and  Mackey-Glass  systems  for  which 
the  fundamental  frequency  is  always  close  to  1/(2 td). 

The  fundamental  frequency  can  be  predicted  from  the  model  through  linear 

stability  analysis  of  the  steady  state.  The  analysis  I  now  present  was  first  performed 

by  my  collaborator,  Dr.  Lucas  Illing.  I  consider  a  small  perturbation  about  the 
steady  state  of  the  form 

P(t)  =  P*  +  6Pext,  (3.69) 

V(t)  =  K  +  6Vext,  (3.70) 

where  SP  and  6V  are  small  amplitudes  and  A  is  the  eigenvalue.  Inserting  these 
solutions  into  Eqs.  3.50  and  3.54,  I  obtain  the  characteristic  equation 

A2  +  (■ —  +  — )  A  +  — -  —  (1  +  baP0)  \e~XTD  =  0,  (3.71) 

\TL  Th  )  TLTh  Tl 

where  I  have  assumed  abP  «  1.  The  time  delay  produces  an  exponential  term 
making  the  characteristic  equation  transcendental.  Since  I  cannot  be  solve  this 
equation  analytically,  I  determine  the  fundamental  frequency  at  the  bifurcation 
point  only.  Assuming  A  =  is  and  separating  the  real  and  imaginary  parts  of  Eq. 
3.71,  I  obtain 

— s2  H - —  (1  +  baPo)  ssin  (std)  —  0,  (3.72) 

tlTh  tl 

( - 1 - )s - —  (1  +  baP0)  s  cos  (std)  =  0.  (3.73) 

\tl  Th  J  tl 

61 


xD  (ns) 


F  igure  3. 17:  Hopf  frequency  as  a  function  of  delay  time  T/>  Circles  are  data  points 
from  the  numerical  model.  Triangles  are  experimental  data  points.  The  dotted  lines 
show  the  function  n / r  where  n  =  1,  2,  3.... 


I  set  x  —  std  and  eliminate  e  from  these  equations  to  obtain 

x 2 - — — f-  td  (  —  +  —  ^  re  tan  (re)  =  0.  (3.74) 

TlTh  \Tl  Th  ) 

This  equation  can  be  solved  numerically  to  find  all  frequencies  that  go  unstable  as 
7  is  increased.  Due  to  the  transcendental  nature  of  the  equation,  there  is  an  infinite 
number  of  such  frequencies.  The  Hopf  frequency  is  the  one  that  goes  unstable  at 
the  lowest  value  of  7.  By  eliminating  773  from  Eqs.  3.73  and  3.73  I  obtain 


th 


1  +  TL/tH 
k  (1  4-  baP0 )  cos  ad 


(3.75) 


which  is  the  threshold  value  of  7  for  a  given  value  of  x.  Given  a  set  of  solutions 


to  Eq.  3.74,  I  use  this  equation  to  find  which  frequency  first  goes  unstable.  The 
result  is  shown  in  by  the  solid  line  in  Fig.  3.18.  As  observed  in  the  experiment, 
the  fundamental  frequency  is  roughly  tl/td  with  discrete  jumps  from  one  branch  to 
another  as  td  increases.  In  the  limit  where  th  — ►  0,  the  same  stability  analysis 


gives  the  fundamental  frequency  the  system  would  display  if  the  high-pass  filter 
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Figure  3.18:  Hopf  frequency  predicted  by  linear  stability  analysis  of  steady  state. 
The  solid  line  corresponds  to  the  active  interferometer  with  bandpass  feedback.  The 
dotted  line  corresponds  to  the  same  system  but  with  the  high-pass  filtering  removed 
from  the  feedback  loop.  The  discrete  frequency  jumps  due  to  the  high-pass  filter  are 
not  displayed  in  well-known  low- pass  systems  such  as  the  Ikeda  and  Mackey- Glass. 

were  removed  from  the  feedback  loop.  In  that  case  Eq.  3.74  becomes 


x  +  —  tan  ( x )  =  0, 
Tl 


(3.76) 


and  Eq.  3.75  becomes 

1th  =  f1  —r  *  - •  (3.77) 

K  (1  T  b \OlFq)  cosx 

The  resulting  frequency,  roughly  1/(2 td),  is  shown  by  the  dashed  line  in  Fig.  3.18. 
This  case  resembles  the  frequencies  observed  in  the  well-known  Ikeda  and  Mackey- 
Glass  systems  that  have  only  low-pass  filtering  in  the  feedback  loop.  The  discrete 
jumps  in  the  fundamental  frequency  can  therefore  be  attributed  to  high-pass  filtering 
in  the  feedback  loop. 

The  origin  of  the  fundamental  frequency  can  be  understood  intuitively  by  con¬ 
sidering  whether  a  wave  circulating  in  the  feedback  loop  will  reinforce  itself.  If  the 
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feedback  is  positive,  a  wave  will  reinforce  itself  if  an  integer  number  of  wavelengths 
equals  the  propagation  length  of  the  loop.  If  the  feedback  gain  is  negative,  as  in 
the  Ikeda  system,  the  propagation  length  must  be  a  half  integer  number  of  wave¬ 
lengths.  Each  of  these  infinite  number  of  wavelengths  represents  a  different  mode  of 
oscillation.  When  the  gain  in  the  loop  is  larger  than  the  losses  of  a  given  mode,  the 
wave  will  grow  each  time  it  travels  around  the  loop.  Assuming  the  gain  bandwidth 
is  not  perfectly  flat,  one  particular  wavelength  will  reach  this  threshold  first  as  the 
gain  is  increased  from  a  low  level.  In  a  system  with  low-pass  feedback,  the  gain  is 
highest  at  low  frequencies  so  the  mode  with  the  lowest  frequency  is  the  fundamental 
oscillation.  If  the  feedback  is  also  high-pass  filtered,  the  gain  may  be  peaked  near  a 
higher  frequency  mode  so  higher  integer  or  half  integer  multiples  of  1/t.d  may  be¬ 
come  the  fundamental  oscillation.  The  Ikeda  and  Mackey-Glass  systems  both  have 
negative,  low-pass  filtered  feedback  so  the  fundamental  frequency  is  always  near 
1/(2 td).  The  active  interferometer  has  positive,  band-pass  filtered  feedback  so  the 
fundamental  frequency  is  roughly  some  integer  multiple  of 1/td  where  the  integer  is 
determined  by  the  shape  of  the  pass  band.  In  the  limit  Th  — »  oo,  the  high-pass  filter 
is  removed  but  a  sign  change  in  the  feedback  is  also  removed  producing  negative, 
low-pass  filtered  feedback  and  consequently  a  fundamental  frequency  of  1/(2td). 

3.4.2  Route  to  Chaos 

Beyond  the  Hopf  bifurcation,  the  system  follows  a  sequence  of  bifurcations  to  chaos. 
The  sequence  does  not  appear  to  follow  any  of  the  standard  routes  such  as  period 
doubling  or  intermittency.  Several  non-standard  routes  to  chaos  have  been  observed 
in  delay  dynamical  systems  including  an  intermittency  route  with  unusual  laminar 
phase  scaling  [102]  and  a  quasi-periodic  route  with  two  period  doublings  [103]. 

Figure  3.19  shows  a  series  of  power  spectra  from  the  experiment  as  7  is  increased. 
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Experimentally  Observed  Behavior 

Feedback  Gain,  7  (mV/mW) 

steady  state 

<5.1±0.5 

period-1 

5.1-11.0  ±1.0 

period-2  to  chaos 

11-17  ±1 

very  broadband  chaos 

>17  ±1 

Table  3.2:  Dynamic  transitions  in  the  opto-electronic  device  without  modulation 
as  the  feedback  gain  is  varied. 

I  measure  the  power  spectral  density  when  the  gain  is  just  below  the  Hopf  bifurcation 
to  obtain  an  estimate  of  the  noise  floor  for  these  and  all  subsequent  frequency  domain 
measurements.  The  resulting  noise  floor  is  approximately  2x1 0 '  3  mW2/MHz2. 
Spectral  features  weaker  than  this  level  are  completely  obscured  in  the  experiment. 

At  feedback  gains  higher  than  the  value  at  the  Hopf  bifurcation,  the  sinusoidal 
periodic  oscillation  begins  to  square  off.  Figure  3.19(a)  shows  a  time  series  and 
power  spectrum  taken  with  7  =  9.4  mV/mW.  The  square  shape  of  the  waveform  is 
most  clearly  indicated  by  the  prominence  of  the  odd  harmonics  in  the  spectrum.  As 
the  feedback  gain  is  increased,  a  small,  broad  peak  appears  at  half  the  fundamental 
frequency  suggesting  a  period-doubling  bifurcation  as  shown  in  Fig.  3.19(b)  where 
7  =  13.2  mV/mW.  The  peak  at  half  the  fundamental  frequency  is  three  orders  of 
magnitude  below  the  fundamental,  indicating  that  the  period-doubled  orbit  only 
deviates  slightly  from  the  original  orbit.  The  weakness  and  broadness  of  the  peak 
at  half  the  fundamental  frequency  coupled  with  the  presence  of  phase  noise  may 
explain  why  the  period  doubling  is  not  apparent  in  the  time  domain  in  Fig.  3.19(b). 
As  the  gain  is  further  increased,  the  broad  background  rises  and  the  tall  peaks  at 
the  fundamental  frequency  and  its  harmonics  weaken.  Figure  3.19(c)  shows  a  time 
series  and  power  spectrum  where  7  =  17.6  mV/mW.  The  power  spectrum  is  quite 
broad  and  the  peaks  have  nearly  dropped  to  the  level  of  the  background.  Table  3.2 
summarizes  the  route  to  chaos.  Although  the  locations  of  bifurcations  are  difficult 
to  determine  precisely  from  the  experimental  system,  the  approximate  range  of  each 
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Figure  3.19:  Experimentally  measured  time  series  and  power  spectra  measured 
from  the  second  output  port  of  interferometer.  Three  steps  on  the  route  to  chaos 
are  shown. 
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type  of  observed  behavior  is  given. 

The  degree  of  complexity  of  the  chaotic  oscillations  is  an  important  quantity  for 
understanding  both  the  dynamics  of  the  system  and  its  usefulness  in  applications. 
For  example,  the  degree  of  security  in  chaos-based  communications  schemes  is  be¬ 
lieved  to  be  related  to  the  complexity  of  the  dynamics  [88]- [90].  Although  techniques 
exist  for  extracting  various  measures  of  complexity  (dimension,  entropy,  Lyapunov 
exponents)  from  experimental  time  series,  they  require  long  (>106  points),  high 
resolution  (12  bit  or  better)  data  sets  with  very  little  noise  [104,  105].  Such  data 
is  not  attainable  from  this  system  due  to  its  high  speed  and  the  strong  influence  of 
phase  noise.  I  have  attempted  to  apply  the  false  nearest  neighbors  [105]  method  for 
estimating  dimension  to  experimental  data  collected  with  a  high-speed  digital  oscil¬ 
loscope.  With  data  sets  of  15,000  points  with  8  bit  resolution,  the  method  failed  to 
converge.  Furthermore,  even  with  a  mathematical  model,  accurate  determination 
of  important  dynamical  quantities  is  difficult  for  any  time-delay  system  exhibiting 
high  dimensional  chaos.  Studies  of  time  delay  systems  with  low-pass  filtered  feed¬ 
back  suggest  that  the  dimension  of  the  chaotic  attractor  is  given  roughly  by  the 
number  harmonics  of  the  fundamental  frequency  that  fit  within  the  low-pass  filter 
bandwidth  in  the  long-delay  limit  [108] .  Loosely  speaking,  increasing  the  time  delay 
puts  more  modes  within  the  bandwidth  of  the  system  increasing  the  complexity  of 
the  dynamics.  Applying  this  relation  to  my  system  gives  an  estimated  maximum 
dimension  of  ~4.5.  However,  it  is  not  clear  that  this  rule  applies  to  this  system. 
Some  attempts  have  been  made  to  estimate  the  Lyapunov  dimension  of  a  delay  sys¬ 
tem  with  band-pass  filtered  feedback  [44],  Interestingly,  the  dimension  was  found 
to  be  nearly  three  times  the  dimension  of  a  low-pass  filtered  system  with  the  same 
bandwidth.  However,  the  value  obtained  (3,700)  is  extremely  high  and  the  method 
used  to  estimate  it  is  generally  believed  to  be  unreliable  for  dimensions  above  10  or 
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so. 

A  similar  route  to  chaos  is  displayed  by  the  model  as  seen  in  Fig.  3.20  where 
time  series  and  power  spectra  are  shown  for  the  values  of  7  used  the  experiment. 
The  power  spectra  are  obtained  by  performing  a  Fast  Fourier  Transform  on  the 
model  timeseries  with  a  610  kHz  frequency  resolution  and  20  GHz  Nyquist  critical 
frequency  [101].  At  7  =  9.4  mV/mW,  the  oscillation,  shown  in  Fig.  3.20(a)  is 
nearly  a  square  wave.  Odd  harmonics  dominate  the  power  spectrum.  At  7  =  13.2 
mV/mW,  the  square  wave  has  undergone  at  least  one  period  doubling;  however,  the 
fundamental  frequency  is  still  several  orders  of  magnitude  larger  than  any  harmonic 
or  subharmonic,  as  shown  in  Fig.  3.20(b).  Also,  a  low  frequency  component  appears 
in  the  time  domain  as  small  differences  in  the  waveform  from  one  period  to  the 
next.  Figure  3.21  shows  the  same  power  spectrum  with  a  hirer  frequency  resolution 
(152.5  kHz).  A  relatively  large  peak  appears  at  1.8  MHz.  This  low  frequency 
oscillation  does  not  appear  to  be  a  transient  or  the  product  of  a  particular  integration 
algorithm.  Since  1.8  MHz  falls  below  the  high-pass  filter  cutoff  frequency  (3  dB 
point),  the  origin  of  this  oscillation  remains  unclear.  The  presence  of  this  low 
frequency  and  its  harmonics  may  contribute  to  the  broadness  of  the  peaks  observed 
in  the  lower  resolution  experimental  power  spectra  at  half  the  fundamental  frequency 
as  for  example  in  Fig.  3.19(b). 

Finally,  at  7  =  17.6  mV/mW,  the  waveform  no  longer  resembles  a  square  wave. 
The  power  spectrum,  shown  in  Fig.  3.20(c),  is  very  broad  with  only  a  few  peaks 
rising  more  than  an  order  of  magnitude  above  the  background.  The  dynamics  of 
the  model  appear  to  resemble  qualitatively  the  experiment  experimentally  observed 
behavior  quite  well. 
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Figure  3.20:  Time  series  and  power  spectral  densities  from  model  at  the  same 
gains  as  experimental  data  shown  in  Fig.  3.19. 
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Figure  3.21:  A  closer  look  at  the  low-frequency  part  of  the  power  spectral  density 
of  the  interferometer  output  from  the  model  with  7  =  13.2  mV/rnW.  The  frequency 
resolution  is  four  times  higher  than  in  Fig.  3.20(b). 

3.5  Dynamics  with  External  Modulation 

In  this  section,  I  investigate  the  effects  of  external  modulation  on  the  active  inter¬ 
ferometer  with  delayed  bandpass  feedback.  First,  I  examine  frequency  locking  due 
to  weak  modulation  near  the  fundamental  frequency.  Frequency  locking  has  been 
reported  in  delay  systems  with  only  low-pass  filtered  feedback  [106,  107]  but  has  not 
previously  been  investigated  in  systems  with  bandpass  feedback.  Frequency  locking 
plays  an  important  role  in  applications  such  as  dynamical  memory  where  the  time- 
delay  system  is  locked  to  an  external  reference  clock  [34],  Next,  I  study  the  route 
to  chaos  under  stronger  external  modulation.  As  shown  in  the  previous  section, 
the  period  doubling  bifurcation  displayed  by  the  active  laser  interferometer  without 
modulation  is  rather  subtle.  The  frequency  component  at  half  the  fundamental  fre¬ 
quency  is  three  orders  of  magnitude  smaller  than  the  component  at  the  fundamental 
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frequency.  In  this  section,  I  show  that  the  addition  of  an  external  modulation  near 
the  fundamental  frequency  significantly  enhances  the  frequency  component  corre¬ 
sponding  to  the  period  doubled  solutions.  This  effect  is  of  practical  importance  to 
the  control  experiments  reported  in  the  next  chapter.  In  control  experiments  gener¬ 
ally,  the  domain  of  control  is  determined  by  varying  the  controller  parameters  and 
observing  an  error  signal  that  measures  the  distance  from  the  desired  orbit.  A  large 
jump  in  the  error  signal  is  interpreted  as  the  edge  of  the  domain.  However,  this 
approach  becomes  problematic  when  applied  to  the  active  interferometer  without 
modulation  because  the  uncontrolled  state  differs  only  slightly  from  the  controlled 
state  so  no  sudden  jump  is  observable.  Since  in  the  modulated  system  the  orbits  are 
more  distinct,  the  error  signal  is  much  easier  to  interpret  and  is  therefore  a  better 
testbed  for  evaluating  the  effectiveness  of  a  controller. 

The  modulation  is  added  by  injecting  a  sinusoidal  voltage  into  the  feedback  loop 
just  before  the  bias-T  through  an  RF  power  combiner  as  shown  in  Fig.  3.3.  The 
effect  of  this  injected  voltage  is  described  by  adding  a  time  dependent  term  to  Eq. 
3.53  to  get 


P(t)  =  crR  [U ( t )  +  Vm  sin  (Qmt)]  +  P0,  (3.78) 


where  Vm  and  Qm/27T  are  the  modulation  amplitude  and  frequency.  Taking  the 
derivative  with  respect  to  time  and  substituting  Eq.  3.52,  the  full  model  becomes 


rLV  =  -V  +  yP^  {1  +  bi  sin  [a  (P^  -  P0)]}  , 


P= - {P-P0)+K 

Th 


V  T  ^lmVm  cos 


(3.79) 

(3.80) 


When  the  modulation  frequency  is  close  to  the  fundamental  frequency,  frequency 
locking  occurs  pulling  the  natural  oscillation  frequency  towards  the  modulation  fre¬ 
quency.  Figure  3.22  shows  the  strongest  peak  in  the  power  spectrum  as  a  function 
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Figure  3.22:  Frequency  locking  due  to  weak  external  modulation  at  a  frequency 
near  the  fundamental  oscillation  frequency.  The  circles  show  experimental  data. 
The  solid  line  shows  results  from  the  mathematical  model.  The  feedback  gain  is 
7  =  6.3  mV/mW. 

of  the  modulation  frequency  from  both  the  experiment  (circles)  and  the  numerically 
integrated  model  (solid  line).  With  the  loop  gain  7  =  6.3  mV/mW  and  the  modu¬ 
lation  amplitude  Vm  =  34  mV,  a  locking  region  of  about  6  MHz  is  apparent.  The 
strong  agreement  between  the  model  and  the  experiment  confirms  that  Eq.  3.80 
correctly  describes  the  effect  of  external  modulation. 

Under  stronger  modulation,  the  system  follows  a  route  to  chaos  much  like  that 
described  in  the  previous  section  where  a  fundamental  frequency  period  doubles  and 
then  broadens  into  chaos  as  7  is  increased.  Figure  3.23  shows  several  experimental 
time  series  and  power  spectra  of  the  optical  power  at  the  second  interferometer  port 
with  Vm  =  225  mV  and  Q,m/2n  =  51.7  MHz.  With  the  feedback  gain  below  7  =  6.7± 
0.4  mV/mW,  the  system  displays  a  periodic  oscillation  at  the  external  modulation 
frequency  (to  within  the  300  kHz  resolution  bandwidth  of  the  spectrum  analyzer) 
as  shown  in  Fig.  3.23(a).  When  the  feedback  gain  is  increased  beyond  7  =  6.7 ±0.4 
mV/mW  a  period  doubling  occurs  producing  a  large  frequency  component  at  half 
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Figure  3.23:  Experimentally  measured  time  series  and  power  spectra  of  the  ouput 
power  at  the  second  port  of  the  active  interferometer  with  external  modulation 
showing  route  to  chaos  as  7  is  increased.  The  loop  gain  7  is  (a)  2.2  mV/mW,  (b) 
8.8  mV/mW,  (c)  11.0  mV/mW,  and  (d)  15.4  mV/mW.  The  period-doubled  orbit  in 
(b)  has  a  much  stronger  component  at  half  the  fundamental  frequency  (~  26  MHz) 
than  in  the  undriven  system  (see  Fig.  3.19(b)). 


Experimentally  Observed  Behavior 

Feedback  Gain,  7  (mV/mW) 

period-1 

<6.7±0.4 

period-2 

6.7-10.0  ±0.5 

period-4  dominated  chaos 

10.0-15.0  ±0.5 

very  broadband  chaos 

>15  ±0.5 

Table  3.3:  Experimentally  observed  dynamic  transitions  in  the  modulated  system 
as  the  feedback  gain  is  varied. 


the  modulation  frequency.  Figure  3.23(b)  shows  the  period  doubled  oscillation  and 
the  large  frequency  component  at  fhn/4vr  when  7  =  8.8  mV/mW.  The  height  of 
the  frequency  component  corresponding  to  the  doubled  period  is  within  an  order 
of  magnitude  of  the  peak  at  the  fundamental  frequency.  This  contrasts  strongly 
with  the  period  doubled  solution  in  the  undriven  system  where  the  period  doubled 
frequency  was  three  orders  of  magnitude  smaller  than  the  fundamental  (see  Fig. 
3.19(b)).  As  the  gain  is  increased  beyond  7  =  10.0  ±  0.5  mV/mW,  the  frequency 
content  around  one  fourth  of  the  fundamental  begins  to  peak  as  seen  in  Fig.  3.23(d) 
where  7  =  11.0  mV/mW.  This  may  indicate  another  period  doubling  transition. 
Finally,  beyond  7  =  15.0  ±  0.5  mV/mW,  the  background  in  the  power  spectrum 
rises  and  the  peaks  fall,  resulting  in  a  broad,  flattened  spectrum  where  the  time  series 
is  extremely  irregular.  The  entire  route  to  chaos  is  summarized  in  Table  3.3.  Using 
the  same  gain  values  as  in  the  experiment,  very  similar  time  series  and  power  spectra 
are  produced  by  the  model,  as  shown  in  Fig.  3.24.  With  7  =  2.2  mV/mW,  as  shown 
in  Fig.  3.24(a),  the  model  displays  a  period  one  orbit.  Around  7  =  6.7  mV/mW, 
a  period  doubling  occurs.  As  in  the  preceeding  section,  the  period-doubled  orbit 
is  subject  to  low-frequency  modulation  as  shown  in  Fig.  3.24(b).  As  the  feedback 
is  further  increased,  the  background  in  the  power  spectrum  broadens  and  a  peak 
appears  at  one  fourth  the  fundamental  frequency,  as  shown  in  Fig  3.24(c).  As  in 
the  experiment,  the  peak  at  one  fourth  the  fundamental  oscillation  is  small  and 
broad.  This  may  indicate  that  the  system  is  already  chaotic  but  that  the  period 
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Figure  3.24:  Time  series  and  power  spectra  from  the  numerically  integrated  model 
of  the  modulated  system  showing  the  route  to  chaos.  As  in  the  experiment,  the 
frequency  content  of  the  period  doubled  orbit  has  a  much  stronger  component  at 
half  the  fundamental  frequency  than  the  similar  orbit  in  the  unmodulated  case  (see 
Fig.  3.20(b)).  The  gain  7  is  (a)  2.2  mV/mW,  (b)  8.8  mV/mW,  (c)  11.0  mV/mW, 
and  (d)  15.4  mV/mW.  75 


four  orbit  is  only  weakly  unstable.  In  Fig.  3.24(d),  where  7  =  2.2  mV/mW,  the 
power  spectrum  is  considerably  flattened  and  the  system  is  clearly  chaotic. 

3.6  Discussion 

In  this  chapter  I  have  presented  a  new  fast  chaotic  time-delayed  optical  system.  As 
shown  in  later  chapters,  the  system  is  suitable  for  fundamental  studies  of  controlling 
fast  dynamics.  It  also  has  many  features  that  make  it  well  suited  for  applications 
such  as  chaos-based  communications. 

The  system  is  flexible  in  the  sense  that  key  characteristics  are  easily  manipulated. 
The  source  of  nonlinearity  in  the  system  is  the  Mach-Zehnder  interferometer.  The 
shape  of  the  nonlinearity  is  easily  modified  by  adjusting  the  mirrors  that  make 
up  the  interfemeter.  In  applications  such  a  synchronization-based  communications, 
where  a  slave  oscillator  must  be  closely  matched  to  a  master,  this  accessibility  is  a 
significant  advantage  over  chaotic  optical  systems  that  are  strongly  affected  by  fixed 
or  non-adjustable  nonlinearities,  such  as  the  material  properties  of  a  semiconductor 
laser  or  erbium-doped  fiber  laser. 

The  timescale  of  the  chaotic  oscillations  is  also  easily  modified.  The  data  and 
analysis  in  Sec.  3.4.1  show  that  the  fundamental  frequency  of  oscillation  can  be 
tuned  by  adjusting  the  feedback  delay  time.  The  total  bandwidth  of  the  chaotic 
fluctuations  is  limited  by  the  passband  of  the  feedback  loop.  For  the  purposes  of 
this  study,  I  choose  to  use  components  that  produce  a  relatively  slow  fundamental 
frequency  of  51  MHz  and  a  pass  band  of  7-250  MHz  so  that  the  system  does  not 
oscillate  too  rapidly  for  a  fast  data  aquisition  device  (a  1  GHz  bandwidth  digitiz¬ 
ing  oscilloscope)  to  obtain  accurate  time  series.  As  shown  in  Fig.  3.17,  by  simply 
reducing  the  feedback  delay  to  8.6  ns,  I  obtain  a  fundamental  oscillation  frequency 
of  103  MHz,  much  faster  than  many  chaotic  opto-electronic  systems  that  have  been 
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reported  in  the  literature  [42,  43,  45].  However,  the  system  could  easily  be  imple¬ 
mented  with  faster  electronic  components  and  a  shorter  feedback  delay  to  obtain 
the  much  higher  speeds  necessary  in  a  real  communication  system.  The  timescale 
of  chaotic  oscillations  in  other  sources  of  fast  optical  chaos,  such  as  the  semicon¬ 
ductor  laser  with  optical  feedback  or  erbium-doped  fiber  laser,  is  determined  by  the 
internal  dynamics  of  a  laser  and,  therefore,  is  not  tuneable. 

The  model  presented  here  may  no  longer  apply  if  the  oscillation  speed  approaches 
the  relaxation  oscillation  frequency.  To  retain  the  model,  the  injection  current  could 
be  set  far  above  threshold  to  achieve  the  highest  possible  relaxation  oscillation 
frequency.  A  value  on  the  order  of  20  GHz  is  obtainable  with  some  diode  lasers  so 
the  model  may  be  scaleable  to  oscillations  at  a  few  gigahertz  [96] . 
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Chapter  4 

M  odified  TDAS  Control  for  Time-Delay 
Systems 


In  the  last  chapter,  I  introduced  a  new  fast  chaotic  time-delay  system.  Fast  time- 
delay  optical  and  electro-optical  chaotic  systems  are  likely  to  be  key  elements  in 
chaos-based  applications  such  as  communications,  dynamic  memory,  and  computa¬ 
tion.  However,  successful  application  of  such  devices  requires  a  fundamental  under¬ 
standing  of  techniques  for  manipulating  the  chaotic  dynamics  such  as  control  and 
synchronization.  Several  researchers  have  demonstrated  techniques  for  synchroniz¬ 
ing  fast  chaotic  optical  systems  [28,  31,  58,  109,  110].  Relatively  little  success  has 
been  demonstrated  with  respect  to  the  related  problem  of  control. 

Much  effort  has  focused  on  controlling  semiconductor  lasers  with  optical  feed¬ 
back,  but  with  limited  success  [38,  54,  111,  112],  The  primary  source  of  difficulty  in 
controlling  such  systems  is  the  high  speed  of  fluctuations.  Rapid  oscillations  present 
two  challenges.  First,  high  speed  sampling  technology  is  necessary  to  obtain  accu¬ 
rate  time  series.  At  present,  the  fastest  digital  sampling  oscilloscopes  can  record 
data  at  rates  of  a  few  billion  samples  per  second.  However,  some  chaotic  optical 
systems  are  known  to  display  oscillations  on  a  time  scale  three  orders  of  magnitude 
shorter.  Second,  as  discussed  in  Ch.  2,  chaos  control  schemes  fail  when  the  control 
loop  latency  is  comparable  to  the  expansion  rate  in  phase  space.  Although  latency 
as  short  as  4.4  nanoseconds  has  been  demonstrated,  this  is  still  clearly  too  long  for 
many  systems  of  interest  [63]. 

Time  delay  auto-synchronization  has  proven  to  be  an  effective  control  technique 
for  chaotic  systems  displaying  oscillations  at  megahertz  frequencies.  In  this  scheme, 
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the  control  signal  is  determined  by  comparing  the  state  of  the  system  at  the  present 
time  with  the  state  of  the  system  at  a  time  in  the  past.  The  state  of  the  system  is 
represented  by  either  one  of  the  variables  of  the  system  or  a  linear  combination  of 
multiple  variables.  The  time  delay  t>  is  chosen  to  be  the  period  7>  of  the  unstable 
periodic  orbit  to  be  stabilized.  Nonetheless,  it  can  fail  if  control  loop  latency  is  larger 
than  10-20%  of  the  period  of  the  orbit.  Extended  time  delay  auto-synchronization 
has  been  shown  to  be  more  robust  to  latency  by  retaining  information  from  an 
infinite  number  of  past  states  in  order  to  determine  the  control  signal.  However, 
implementing  ETDAS  can  be  very  challenging. 

In  this  chapter,  I  introduce  a  modified  approach  to  TDAS  control  applicable  to 
time  delay  systems  that  allows  for  latency  equal  to  the  time  delay.  This  technique 
should  be  applicable  to  many  of  the  electro-optic  time  delay  systems  currently  be¬ 
ing  studied  and  may  also  be  useful  for  other  systems  described  by  delay  differential 
equations.  I  apply  the  modified  TDAS  control  approach  to  the  active  laser  inter¬ 
ferometer  with  bandpass  feedback  described  in  the  previous  chapter.  In  Sec.  4.1,  I 
describe  this  new  control  approach.  In  Sec.  4.2,  I  present  an  experimental  demon¬ 
stration  of  the  modified  TDAS  approach  by  applying  it  to  the  active  interferometer 
with  bandpass  feedback  and  external  modulation.  The  success  of  control  is  quan¬ 
tified  by  measuring  the  domain  of  control.  In  Sec.  4.3,  I  develop  a  mathematical 
model  of  the  interferometer  with  control.  The  model  shows  behavior  in  qualitative 
agreement  with  the  experiment. 

4.1  M  odified  TDAS  Control 

Several  optical  time-delay  systems  have  been  designed  expressly  to  produce  sub¬ 
nanosecond  chaotic  fluctuations  [27,  44,  57].  Such  fast  fluctuations  would  be  difficult 
to  control  using  techniques  that  cannot  compensate  for  a  significant  amount  of 
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control  loop  latency.  In  this  section,  I  introduce  a  modified  form  of  the  TDAS  control 
method  for  fast  chaotic  time-delay  systems  that  permits  an  amount  of  latency  equal 
to  the  feedback  delay  time  of  the  system. 

To  illustrate  this  control  approach,  I  will  apply  it  the  model  equation  introduced 
by  Ikeda 

x(t)  =  -x(t)'*jr'yf[x(t-TD)\,  (4.1) 


where  7  is  the  strength  of  the  nonlinearity  and  Td  is  the  delay  time.  The  function 
f  (x)  is  assumed  to  have  at  least  one  maximum.  This  equation  was  originally  derived 
as  a  model  of  a  ring  cavity  containing  a  nonlinear  dielectric  medium  in  which  case 
f  [x\  =  7r  sin(x)  [76].  Recently,  numerous  electro-optic  devices  have  been  devised 
that  are  described  by  this  type  of  equation  with  various  nonlinear  functions  [43], 
[110]-[35] .  The  Mackey-Glass  Equation,  a  common  model  for  studying  time  delay 
systems,  is  also  of  this  form  with  /  [x]  =  x(l  +  x10)-1[41]. 

In  all  these  cases,  Eq.  4.1  may  be  viewed  as  describing  a  low  pass  filter  with 
unit  response  time  and  input  u(t)  =  7/  [x  (t  —  td)\  as  shown  schematically  in  Fig. 
4.1.  The  input  term  u(t)  retains  a  “memory”  of  the  system  state  a  time  Td  in  the 
past.  This  memory  is  reflected  in  the  autocorrelation  function  of  f  [x  (£)]  defined  as 


Cauto  \J~D) 


f-oo  f[x(t)\f[x(t-TD  )}dt 


(4.2) 


f-oo  f[x(t)]2dt 

which  gives  the  correlation  between  a  time  series  of  f  [x  (£)]  and  the  same  time  series 
with  a  delay  Td  [104],  For  example,  Fig.  4.2  shows  the  autocorrelation  function 
of  the  input  term  of  Eq.  4.1  where  f  [x\  =  7rsin(x)  with  7  =  2.1  and  Td  =  20. 
With  these  parameter  values,  the  dimension  of  the  system  has  been  estimated  to 
be  approximately  42  [40].  Nonetheless,  a  series  of  peaks  in  the  autocorrelation 
function  appear  at  frequencies  near  utd  where  n  =  1,2,3,...  .  This  suggests  that 
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Figure  4.1:  Block  diagram  showing  the  general  form  of  the  class  of  time-delay 
systems  to  which  the  modified  TDAS  method  may  be  applied.  The  basic  elements 
are  a  low  pass  filter,  a  nonlinear  element  (f [x] ) ,  and  a  delayed  feedback  loop.  Well 
known  models  such  as  the  Ikeda  and  Mackey-Glass  equations  are  of  this  form. 

the  f  [x(t  —  Tjj)]  may  contain  sufficient  information  to  predict  f[x(t)]  and  perturb 
the  system  accordingly  in  order  to  effect  control. 

Based  on  this  observation,  I  propose  to  control  time  delay  systems  by  modifying 
the  TDAS  approach  in  two  ways.  First,  the  control  loop  latency  is  intentionally  set 
to  To-  Second,  the  feedback  control  signal  is  the  difference  between  the  values  of 
the  nonlinear  function  /  [x]  at  times  separated  by  the  desired  orbits  period  rather 
than  the  difference  in  the  values  of  x  itself.  Applying  this  control  method  to  Eq. 
4.1  gives 

x(t)  =  -x(t)  +  7 /  [x  (t  -  TD)\  +  7c  {/  [x  (t-TD)\-  /  [x  ( t-To  -  TP )]}  , 

(4.3) 

where  yc  is  the  control  gain  and  7>  is  the  period  of  the  desired  orbit.  Figure  4.3 
shows  a  block  diagram  of  the  system  plus  controller. 

The  original  TDAS  scheme  assumes  nearly  instantaneous  feedback  from  the 
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delay 


Figure  4.2:  Autocorrelation  function  of  the  nonlinear  term  of  the  Ikeda  model  of 
a  passive  nonlinear  resonator  in  the  regime  of  high  dimensional  chaos.  The  peak 
at  td  suggests  that  a  delayed  version  of  the  nonlinear  term  may  contain  sufficient 
information  about  the  current  state  of  the  system  to  successfully  control  it. 

controller,  whereas  it  has  been  shown  that  latency  on  the  order  of  20%  of  the 
period  of  the  orbit  is  usually  enough  to  foil  TDAS  control  [49,  50].  The  fundamental 
oscillation  displayed  by  this  class  of  systems  typically  has  a  period  of  Tp>  or  2 Tp  so 
a  TDAS  controller  is  likely  to  fail  if  subject  to  a  latency  of  Tp.  The  modified  TDAS 
scheme  allows  a  period  of  time  Tp  in  order  to  compute  and  apply  the  appropriate 
control  signal  and  therefore  may  allow  control  of  faster  systems. 

To  test  the  proposed  control  scheme  experimentally,  I  apply  it  to  the  time  delay 
system  described  in  the  previous  chapter,  the  active  interferometer  with  bandpass 
feedback  and  external  modulation.  The  control  signal  is  generated  using  the  second 
output  of  the  interferometer.  A  photodetector  at  the  second  port  of  the  interferom¬ 
eter  converts  the  emitted  power  into  a  voltage.  The  voltage  signal  is  split.  One  half 
is  delayed  by  a  time  Tp  and  then  subtracted  from  the  other.  The  resulting  signal 
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Figure  4.3:  Block  diagram  showing  the  modified  TDAS  control  scheme.  The 
feedback  control  signal  is  the  difference  between  the  values  of  the  nonlinear  function 
/  [re]  at  times  separated  by  the  period  of  the  desired  orbit.  The  control  loop  latency 
is  set  equal  to  td 

is  amplified  and  combined  with  the  feedback  loop  of  the  interferometer  before  the 
bias-T.  Figure  4.4  shows  a  schematic  diagram  of  the  entire  system  including  the 
controller. 

4.2  Experimental  Demonstration  of  M  odified  TDAS 
Control 


In  this  section,  I  describe  the  experimental  demonstration  of  the  modified  TDAS 
control  scheme.  The  experimental  system  is  the  active  interferometer  with  delayed 
bandpass  feedback  and  external  modulation.  I  describe  how  I  implemented  the  con¬ 
troller  using  radio-frequency  electronic  components.  I  present  methods  for  properly 
setting  the  necessary  time  delays.  The  effectiveness  of  the  controller  is  quantified 
by  estimating  the  domain  of  control  of  a  periodic  orbit  with  a  frequency  close  to 
that  of  the  fundamental  oscillation  frequency. 
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Figure  4.4:  Schematic  diagram  of  active  interferometer  and  controller.  A  pho¬ 
todetector  at  the  second  port  of  the  interferometer  converts  the  emitted  power  into 
a  voltage.  The  voltage  signal  is  split.  One  half  is  delayed  by  a  time  7>  and  then 
subtracted  from  the  other.  The  resulting  signal  is  amplified  and  combined  with 
the  feedback  loop  of  the  interferometer  before  the  bias-T.  The  propagation  time 
through  the  controller  is  set  equal  to  the  feedback  delay  time,  tq. 
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4.2.1  Modified  TDAS  Controller  Implementation 


PF 


PF 


F  igure  4.5:  Schematic  diagram  of  the  modified  TDAS  controller.  The  components 
labelled  A-I  are  as  follows:  A  -  Hamamatsu  silicon  photodiode  S4751;  B  -  Coax¬ 
ial  cable  RG  58/U  and  variable  delay;  C  -  MiniCircuits  power  splitter /combiner 
ZFSC-2-1W-75;  D  -  Coaxial  cable  RG  58/U  and  variable  delay;  E  -  M/ACOM  180 
hybrid  junction  H-9;  F  -  MiniCircuits  directional  coupler  ZFDC-20-4;  G  -  MiniCir¬ 
cuits  amplifier  ZFL-1000GH;  H  -  MiniCircuits  directional  coupler  ZFDC-10-1;  I  - 
MiniCircuits  power  splitter/combiner  ZFSC-2-1W-75. 


To  implement  the  modified  TDAS  scheme,  I  use  RF  electronic  components  sim¬ 
ilar  to  those  used  in  the  feedback  loop  of  the  active  interferometer  itself.  Figure  4.5 
shows  the  experimental  implemenation  of  the  modified  TDAS  controller  in  detail. 
The  intensity  at  the  output  of  the  second  port  of  the  interferometer  falls  on  a  Si 
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photodiode  (A  in  Fig.  4.5).  Only  23%  of  the  power  is  detected  because  the  beam 
width  is  larger  than  the  active  region  of  the  detector  for  reasons  discussed  in  Sec. 
3.1.  The  photodiode  produces  a  current  proportional  to  the  incident  power  that  is 
converted  to  a  voltage  through  a  50  O  resistance. 

This  voltage  signal  propagates  through  a  length  of  coaxial  cable  to  an  RF  power 
splitter  (C  in  Fig.  4.5,  5-600  MHz  bandwidth).  The  length  of  the  cable  is  chosen 
so  that  the  total  propagation  time  through  the  entire  control  loop  is  td.  (The 
procedure  for  accurately  determining  the  proper  length  is  described  below.)  The 
voltage  at  each  of  the  two  outputs  of  the  power  splitter  is  l/\/2  times  the  input 
voltage.  The  voltage  signal  from  port  1  of  the  power  splitter  travels  directly  to  the 
180°  hybrid  junction  (M/ACOM  H-9,  2-2000  MHz  bandwidth,  labelled  E  in  4.5). 
The  voltage  signal  from  port  2  of  the  splitter  travels  down  a  length  of  coaxial  cable 
before  reaching  the  hybrid  junction.  The  length  of  the  cable  (D  in  Fig.  4.5)  is 
chosen  so  that  the  voltage  signal  from  port  2  is  delayed  by  a  period  of  time  equal 
to  Tp  relative  to  the  voltage  signal  from  port  1. 

The  two  voltage  signals  from  the  splitter  are  recombined  by  the  180°  hybrid 
junction  which  performs  the  subtraction  operation.  This  passive  device  has  two 
input  ports,  labelled  +  and  —  in  Fig.  4.5,  and  two  output  ports  labelled  A  and  E. 
When  RF  voltages  V+  and  VI  are  applied  to  the  +  and  —  ports  respectively,  the 
voltages  at  the  A  and  E  ports  are 


^(V-vq, 

(4.4) 

7!(Vl  +  v-h 

(4.5) 

The  output  of  port  1  of  the  power  splitter  enters  the  +  port.  The  signal  from  port 
2  of  the  power  splitter  enters  the  —  port  of  the  hybrid  after  propagating  through 
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the  cable.  The  E  port  of  the  hybrid  is  terminated  through  50  and  the  power 
dissipated  here  is  not  used.  The  A  port  supplies  the  signal  that  travels  through  the 
remainder  of  the  control  loop. 

A  small  fraction  (1%)  of  the  power  output  from  the  A  port  of  the  hybrid  is  sam¬ 
pled  by  a  directional  coupler  (MiniCircuits  ZFDC-20-4,  1-1000  MHz  bandwidth), 
amplified  by  a  low  noise  amplifier  (MiniCircuits  ZFL1000LN,  0.1-1000  MHz  band¬ 
width),  and  recorded  by  a  high-speed  digital  oscilloscope  (Tektronix  TDS680B,  1 
GHz  bandwidth).  This  signal  provides  a  measure  of  the  error  between  the  state  of 
the  system  and  the  desired  state. 

The  remainder  of  the  signal  power  from  the  hybrid  continues  in  the  control 
loop  passing  through  a  220  pF  capacitor,  a  variable  gain  amplifier  (MiniCircuits 
ZFL-1000GH,  10-1000  MHz  bandwidth),  and  a  470  pF  capacitor.  The  capacitors 
are  chosen  to  match  those  in  the  interferometer  feedback  loop  (see  Fig.  3.3)  so 
that  the  controlled  signal  undergoes  the  same  band  pass  filtering  as  the  feedback 
signal.  The  resulting  control  signal  is  combined  with  the  external  modulation  signal 
through  a  directional  coupler  (MiniCircuits  ZFDC-10-1,  1-500  MHz  bandwidth). 
The  control  signal  is  directed  into  what  is  normally  the  “output”  port  while  the 
modulation  signal  enters  through  the  coupling  port.  The  combined  signal  exits 
through  the  nominal  “input”  port.  This  configuration  is  chosen  in  order  to  isolate 
the  components  that  measure  the  error  signal  from  the  modulation  source.  Finally, 
the  combined  modulation  and  control  signal  reaches  a  power  combiner  that  sums  it 
with  the  interferometer’s  feedback  signal.  The  resulting  signal  is  applied  to  the  RF 
input  of  the  bias-T  at  the  laser  mount. 

Both  sections  of  coaxial  cable  used  to  set  delay  times  contain  adjustable  length 
coaxial  lines  (General  Radio  Co.  874-LK10L).  These  devices,  somewhat  resembling 
a  trombone  slide,  allow  fine  tuning  of  the  cable  delay  over  a  range  of  310  ps.  The 
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delay  line  (B  in  Fig.  4.5)  located  immediately  after  the  photodetector  sets  the 
control  loop  latency.  Its  length  is  chosen  to  provide  a  total  latency  19.1  nanoseconds, 
he.,  equal  to  To-  The  second  delay  line  (D  in  Fig.  4.5),  which  sets  the  orbit  period, 
is  adjusted  to  minimize  the  control  error  signal  inside  the  domain  of  control.  A 
control  delay  of  19.33±0.05  ns  appears  to  minimize  the  error  signal  throughtout  the 
experimentally  observed  domain  of  control. 

4.2.2  Procedure  for  Measuring  Delays 

In  order  to  correctly  implement  the  modified  TDAS  control  scheme,  two  path  lengths 
must  be  carefully  determined  in  order  to  obtain  the  appropriate  time  delays  and 
To-  Different  techniques  are  appropriate  for  determining  7>  and  To  since  the  former 
is  due  entirely  to  propagation  through  coaxial  cable  while  the  latter  includes  phase 
shifts  due  to  amplifiers  and  filters,  and  propagation  through  both  free  space  and 
coaxial  cable.  Therefore  I  describe  each  method  separately. 

The  path  that  produces  the  time  delay  Tp  is  the  length  of  cable  between  the 
power  splitter  and  the  -  input  port  of  the  hybrid  minus  the  length  between  the 
splitter  and  the  +  input  port.  The  delay  associated  with  this  path  is  determined  by 
applying  a  sinusoidal  voltage  to  the  input  of  the  splitter  and  observing  the  signal  at 
the  A  port  of  the  hybrid.  The  transfer  function  through  the  splitter-delay-hybrid 
combination  is  derived  by  assuming  an  input  of  Vjn  =  V§ewt  and  a  response  in  the 
form  Vout  =  V/\eluJt  where  Va  is  an  amplitude  to  be  determined.  Inserting  these 
functions  into  Eq.  4.4,  including  a  factor  of  l/y7^  due  to  the  splitter,  and  assuming 


the  delay  of  one  path  relative  to  the  other  is  r  gives 


Va 

Va 

V0 


\  (r0e“  -  Vie1""-'’) , 


(4.6) 

(4.7) 


From  Eq.  4.7,  it  can  be  seen  that  the  transfer  function  has  notches  at  all  integer 
multiples  of  1/r.  The  location  of  the  first  notch  is  determined  by  analyzing  the 
output  of  the  A  port  with  a  spectrum  analyzer  with  resolution  bandwidth  of  30kHz. 
The  resulting  uncertainty  in  the  delay  time  is  on  the  order  of  0.01  ns. 

The  path  that  produces  control  loop  latency  of  Td  contains  the  entire  controller 
including  propagation  through  the  interferometer,  splitter,  hybrid,  cables,  capaci¬ 
tors,  and  amplifier.  Many  of  these  components  add  phase  delay  in  addition  to  prop¬ 
agation  delay.  Therefore,  a  direct  measurement  of  the  propagation  distance  only 
gives  an  order-of-magnitude  estimate  of  the  delay  time.  To  get  a  more  accurate 
measure  of  the  delay,  I  take  advantage  of  the  fact  that  the  active  interferometer 
without  external  modulation  oscillates  at  a  fundamental  frequency  very  close  to 
1/td.  By  blocking  the  detector  at  the  first  port  of  the  interferometer  and  removing 
the  delay  line  from  between  the  splitter  and  hybrid,  I  create  a  new  feedback  loop.  I 
increase  the  gain  until  oscillations  appear.  The  signal  from  the  directional  coupler 
(labelled  F  in  Fig.  4.5)  is  measured  by  a  spectrum  analyzer  to  determine  the  fre¬ 
quency  of  oscillation.  I  estimate  the  error  in  the  delay  determined  in  this  manner 
to  be  0.2  ns  or  less. 


4.2.3  Observation  of  Control 

I  characterize  the  effect  of  the  controller  on  the  experimental  system  in  two  ways. 
First,  I  examine  the  changes  in  the  power  spectrum  of  the  interferometer  output  as 
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(a) 


Figure  4.6:  Typical  experimental  power  spectra  recorded  when  the  feedback  gain 
7  =  12.4.  (a)  The  control  gain  7C  =  0.08  is  small  enough  that  the  spectrum  is 
indistinguishable  from  that  of  the  system  with  no  control,  (b)  The  control  gain 
7C  =  5.5,  well  within  the  domain  of  control.  The  large  peaks  at  half-integer  mul¬ 
tiples  of  the  fundamental  frequency  are  suppressed  more  than  20  dB  below  their 
uncontrolled  amplitude  and  the  broad  background  falls  several  dB  below  its  un¬ 
controlled  level,  (c)  The  control  gain  yc  =  9.4,  new  sideband  frequencies  appear 
around  the  fundamental  frequency  and  its  harmonics  indicating  control  is  no  longer 
effective. 
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the  control  gain  is  varied.  Then  I  estimate  the  domain  of  control  by  examining  the 
size  of  the  error  signal. 

To  observe  the  influence  of  the  controller  on  the  system  in  the  frequency  domain, 
the  signal  from  the  directional  coupler  in  the  feedback  loop  (labeled  B  in  Fig.  3.3)  is 
amplified  and  sent  to  a  spectrum  analyzer.  Figure  4.6  shows  a  set  of  typical  power 
spectra  recorded  when  the  feedback  gain  7  =  12.4  mV/mW.  In  Fig.  4.6(a),  the 
control  gain  7C  =  0.08  mV/mW,  a  small  enough  value  that  the  spectrum  is  indis¬ 
tinguishable  from  that  of  the  system  with  no  control.  Dominant  peaks,  located  at 
integer  and  half-integer  multiples  of  the  fundamental  frequency,  rise  out  of  a  broad 
background  indicating  chaotic  oscillation.  In  Fig.  4.6(b),  7C  =  5.5  mV/mW.  The 
large  peaks  at  half-integer  multiples  of  the  fundamental  frequency  are  suppressed 
more  than  20  dB  below  their  uncontrolled  amplitude.  Also,  the  broad  background 
falls  several  dB  below  its  uncontrolled  level.  The  system  at  this  point  is  well  within 
the  domain  of  control.  In  Fig.  4.6(c),  where  7 c  =  9.4  mV/mW,  new  sideband 
frequencies  appear  around  the  fundamental  frequency  and  its  harmonics  indicating 
control  is  no  longer  effective.  The  appearance  of  sidebands  as  the  control  gain  is 
increased  to  a  value  beyond  the  domain  of  control  has  been  observed  in  an  experi¬ 
mental  study  of  the  original  TDAS  method  [49]. 

A  broader  view  of  the  effect  of  the  controller  is  obtained  by  measuring  the  domain 
of  control.  To  do  this,  I  use  the  signal  from  the  directional  coupler  labelled  F  in 
Fig.  4.5  as  a  measure  of  the  control  error  signal.  When  the  feedback  gain  is  set 
to  a  low  value  where  the  system  oscillates  at  the  fundamental  frequency,  the  root 
mean  square  value  of  the  error  signal  is  3.5  ±1  millivolts.  This  signal  is  due  to  noise 
from  the  laser  and  amplifiers  and  to  imperfections  in  the  controller.  I  consider  this 
level  a  good  estimate  of  the  error  signal  I  can  expect  when  the  system  is  within 
the  domain  of  control  and  choose  the  cutoff  value  to  be  just  slightly  larger  at  4.2 
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Figure  4.7:  Experimentally  measured  average  error  signal  at  three  different  gain 
values.  The  dashed  line  shows  the  cutoff  of  4.2  mV  chosen  to  define  the  domain  of 
control. 

mV.  I  vary  the  feedback  gain  and  control  gain  over  a  wide  range,  recording  the  root 
mean  square  value  of  the  error  signal  at  several  points.  Those  points  where  the 
error  signal  average  is  below  the  cutoff  are  considered  to  be  within  the  domain  of 
control. 

Data  collection  is  automated  using  a  personal  computer  and  digital  oscilloscope. 
The  computer  and  oscilloscope  communicate  through  a  GPIB  interface.  The  com¬ 
puter  contains  two  analog  outputs  that  are  used  to  set  7  and  jc.  After  setting  the 
gain  values,  the  computer  directs  the  digital  oscilloscope  to  record  a  time  series  of 
2,500  points  at  a  sample  rate  of  500  MS/s  and  calculate  the  root  mean  square  value. 
This  value  is  relayed  to  the  computer  for  storage  and  the  gains  are  adjusted  to  the 
next  value. 

Three  data  sets  collected  in  this  manner  are  shown  in  Fig.  4.7.  Each  set  is 
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Figure  4.8:  Experimentally  measured  domain  of  control.  At  each  point,  the  root 
mean  square  of  the  error  signal  was  compared  with  a  cutoff  value.  The  black  dots 
mark  points  below  the  cutoff. 

recorded  with  a  different  feedback  gain.  First,  the  circles  show  the  error  signal  when 
7  =  8.5  mV/mW.  At  this  gain  the  uncontrolled  system  is  just  above  the  transition 
from  oscillation  at  the  fundamental  frequency  to  oscillation  at  half  the  fundamental 
frequency.  Over  a  large  range  of  yc,  the  error  signal  is  below  the  cutoff  value  marked 
by  the  dashed  line.  The  triangles  show  the  average  error  at  a  higher  gain  where 
7  =  12.0  mV/mW.  At  this  feedback  gain,  the  uncontrolled  system  displays  oscilla¬ 
tions  with  a  strong  frequency  component  at  half  the  fundamental  frequency  and  a 
broadening  background  including  a  peak  at  one  fourth  the  fundamental  frequency. 
Only  in  a  relatively  small  part  of  the  range  of  7C  does  the  error  signal  remain  below 
the  cutoff.  Finally,  at  a  high  gain  value  of  7C  =  15.7  mV/mW,  where  the  oscillations 
of  the  uncontrolled  system  are  broadband,  the  squares  show  that,  although  the  error 
signal  does  display  a  minimum,  the  error  signal  never  falls  below  the  cutoff.  The 
system  is  considered  to  be  outside  the  domain  of  control  at  these  points. 
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Chaos  Control  M  ethod 

Orbit  Period,  rP 

Control  Loop  Latency,  re 

TDAS  [49,  72] 

99  ns 

<30  ns 

ETDAS  [49,  72] 

99  ns 

<120  ns 

Pulsewidth  Modulation  [63] 

52  ns 

4.4  ns 

Modihed  TDAS 

19.33  ns 

19.1  ns 

Table  4. 1:  Shortest  period  orbits  controlled  experimentally  in  chaotic  systems.  The 
modihed  TDAS  method  described  in  this  chapter  is  used  to  stabilize  the  shortest 
timescale  dynamics  reported  to  date.  Note  that  with  this  method  the  latency  need 
not  be  small  compared  to  the  period  of  the  orbit. 

Having  explained  the  process  by  which  I  quantify  the  effect  of  the  controller,  I 
now  present  my  experimental  estimation  of  the  control  domain.  Figure  4.8  shows 
all  points  within  the  range  of  7  and  yc  at  which  the  error  signal  was  observed  to 
be  below  the  cutoff.  The  estimated  domain  has  a  triangular  shape  such  that  the 
range  of  effective  control  gains  decreases  approximately  linearly  as  the  feedback 
gain  is  increased.  As  might  be  expected,  the  domain  is  very  wide  at  feedback  gain 
values  where  the  uncontrolled  system  is  either  periodic  or  only  slightly  chaotic. 
As  the  feedback  gain  increases  and  the  system  becomes  more  chaotic,  the  domain 
is  thinner  but  still  significant.  Control  apparently  fails  completely  for  feedback 
gain  values  where  the  uncontrolled  system  shows  a  broad  and  nearly  flat  spectrum 
(roughly  7  =  14.5  mV/mW).  A  few  “holes”  appear  at  the  edges  of  the  domain 
where  the  error  signal  is  just  above  the  cutoff.  Such  behavior  is  typical  of  a  noisy 
system  near  a  stability  threshold. 

The  experimentally  observed  domain  of  control  shown  in  Fig.  4.8  extends  from 
the  period-2  region  into  the  chaotic  regime  where  the  period-2  frequency  is  still  quite 
large  noted  in  Table  3.3.  In  the  chaotic  region,  the  modihed  TDAS  control  method 
successfully  stabilizes  a  periodic  orbit  with  a  period  of  19.33  ns.  In  Table  4.1,  this 
result  is  compared  with  the  three  shortest  period  experimentally  stabilized  orbits 
of  chaotic  systems  reported  in  the  literature.  The  period  of  the  orbit  stabilized 
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here  is  less  than  half  the  period  of  the  next  fastest  orbit.  Also  given  in  Table  4.1 
are  the  experimentally  successful  control  loop  latencies  for  each  method.  Note  that 
the  latency  of  the  modified  TDAS  controller  is  nearly  equal  to  the  stabilized  orbit 
period.  Only  one  other  chaos  control  method,  ETDAS,  has  demonstrated  a  similar 
tolerance  to  latency.  Experimental  implementation  of  ETDAS  is  likely  to  be  more 
challenging  than  the  modified  TDAS  method  in  many  cases  because  it  can  require  a 
predistortion  stage  in  the  controller  [49].  Thus,  the  modified  TDAS  method  provides 
a  balance  of  design  simplicity  and  tolerance  to  latency. 

4.3  Mathematical  Model  of  Interferometer  with 
Controller 

To  gain  more  insight  into  the  experimental  results  presented  above,  I  develop  a 
mathematical  model  of  the  controlled  system.  Based  on  the  detailed  description 
of  the  controller  implementation  in  the  last  subsection,  the  interferometer  model 
developed  in  the  previous  chapter  is  now  easily  augmented  to  include  the  modified 
TDAS  controller.  In  this  subsection,  I  add  control  to  the  model  and  then  use  it  to 
predict  the  domain  of  control. 

According  to  Eq.  3.43,  the  power  P2  emitted  by  the  second  port  of  the  interfer¬ 
ometer  is  given  by 

P2(t)  =  *fP(t)  {1  -  b sin  [a  (P(t)  -  P0)]}  ,  (4.8) 

where  P(t)  is  the  power  emitted  by  the  laser  and  b,  a,  Po,  and  (f  are  as  given  in 
Sec.  3.2.3.  A  fraction  Qc  of  this  power  falls  on  the  detector,  which  converts  it  to 
current  with  a  sensitivity  of  p  =  0.44  mA/mW.  A  resistor  r  =  50  D  converts  the 
current  into  the  voltage 

Vdet(t)  =  ^y^P(t)  {1  -  b  sin  [a  (P(t)  -  P0)]}  .  (4.9) 
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This  voltage  is  delayed  and  split  into  two  separate  signals.  One  of  these  is  delayed 
again.  Then  the  two  are  recombined  by  the  hybrid  junction  according  to  Eq.  4.4 
producing  a  voltage 

Vhybrid(t)  =  {P(t)  [1  -  b  sin  [a  (P(t)  -  P0)]]  -  Pt?  [1  -  &  sin  [a  (PTp  -  P0)]]}  , 

4  (4-10) 

where  I  use  the  notation  PT  =  P(t  —  r)  and  I  have  ignored  for  the  moment  all 

other  propagation  delays.  A  directional  coupler  splits  off  1%  of  the  power  from  this 

signal  reducing  the  voltage  by  a  factor  of  0.99.  The  capacitors  and  amplifier  in  the 

control  loop  are  chosen  to  match  the  components  in  the  interferometer  feedback 

loop.  This  is  done  so  that  the  control  signal  receives  the  same  bandpass  filtering  as 

the  feedback  signal.  The  entire  system  is  then  described  by  the  equations 

tlV  =  -V  +  7 Pro  {1  +  b  sin  [a  (PTD  -  P0)]|  +  7 Cc(t),  (4.11) 


P 


1 

Th 


(P  Pq)  +  K 


V  T  ^mVm  COS  (flmt)  , 


(4.12) 


c(t)  -  { {1  -  bsm  [a  (Pn  -  F0)]}  -  Pn+n  {1  -  l.sm  \o  (7VIJ+rp  -  i=o)]}}  • 

(4.13) 

The  quantity  c(t)  will  be  used  below  as  a  measure  of  the  effectiveness  of  the  con¬ 
troller.  The  control  gain  7C  is  given  by 


7c 


P'^CcCfC 

4 


(4.14) 


I  use  this  model  to  obtain  a  theoretical  estimate  of  the  domain  of  control.  For 
the  original  TDAS  scheme,  tools  have  been  developed  to  locate  or  approximate  the 
domain  of  control  that  don’t  require  integration  of  the  full  set  of  equations  governing 
the  system  and  controller.  However,  these  techniques  have  not  been  generalized  to 
time-delay  dynamical  systems  or  to  include  latency.  Therefore,  I  resort  to  direct 
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Figure  4.9:  Time  series  of  power  emitted  from  created  by  numerically  integrating 
model  of  inteferometer  and  controller.  The  feedback  gain  is  7  =  14.0  mV/mW.  (a) 
Uncontrolled  irregular  oscillation  with  yc  =  0  mV/mW.  (b)  Stabilized  oscillation 
near  the  fundamental  frequency  with  7C  =  4.5  mV/mW. 
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numerical  integration  of  Eqs.  4.11-4.13  in  order  to  locate  the  domain  of  control.  I 
use  the  fixed  step-size  an  Adam-Bashforth-Moulton  predictor-corrector  algorithm 
[101]. 

There  are  some  difficulties  with  this  approach.  First,  the  initial  conditions  of 
the  system  consist  of  two  functions  (one  each  for  P  and  V)  defined  on  the  interval 
[— r,  0] .  In  a  phase  space  of  infinite  dimensions,  there  is  no  obvious  method  for 
sampling  a  large  number  of  initial  conditions  in  order  to  determine  where  a  typical 
trajectory  will  evolve.  Moreover,  it  has  been  observed  that  the  basin  of  attraction 
of  a  stabilized  UPO  shrinks  as  the  boundary  of  the  domain  of  control  is  approached. 
Since  only  a  finite  number  of  inital  conditions  can  be  investigated,  this  approach 
necessarily  underestimates  the  true  domain  of  control  to  some  extent.  For  points 
well  within  the  domain  of  control,  I  choose  a  sine  wave  with  the  frequency  of  the 
orbit  to  be  controlled  as  the  initial  condition  of  the  system  with  the  intention  of 
starting  the  system  close  to  the  desired  orbit  in  phase  space.  For  points  closer  to 
the  boundary  of  the  domain  of  control,  I  use  a  time  series  of  the  controlled  orbit 
near  the  center  of  the  domain  of  control  as  the  initial  condition. 

A  second  difficulty  is  that  generally  the  period  of  the  unstable  orbit  Tp  is  not 
commensurate  with  the  feedback  delay,  Tp.  Since  the  integration  step  size  is  chosen 
such  that  an  integer  number  of  steps  fit  within  Tp,  the  same  is  not  be  true  for  Tp. 
The  period  Tp  can  only  be  specified  to  within  the  accuracy  of  the  integration  step 
size.  As  a  result  of  this  error  in  Tp,  the  error  signal  never  goes  to  zero.  In  spite 
of  these  issues,  I  proceed  to  use  numerical  integration  to  estimate  the  domain  of 
control  for  though  it  may  be  a  crude  tool,  it  happens  to  be  the  only  one  available 
at  this  time. 

Typical  results  from  the  integrated  model  are  shown  in  Fig.  4.9.  The  power 
emitted  by  the  laser  is  plotted  as  a  function  of  time  with  and  without  control. 
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Figure  4.10:  Root  mean  square  of  error  signal  from  model.  Domain  of  control 
is  estimated  by  recording  control  gain  values  that  produce  an  error  signal  below  a 
cutoff  value.  The  dashed  horizontal  line  shows  the  cutoff  at  0.01  mW. 

With  7  =  14.0  mV/mW,  the  uncontrolled  system  (see  Fig.  4.9(a))  displays  chaotic 
oscillations  (see  the  preceding  chapter  for  more  description  of  the  uncontrolled  dy¬ 
namics).  When  the  control  is  applied  (see  Fig.  4.9(b)  where  7C  =  4.5)  a  periodic 
oscillation  appears. 

I  adjust  the  control  delay  to  minimize  the  root  mean  square  value  of  the  error 
signal.  In  general,  the  period  of  an  unstable  orbit  in  a  dynamical  system  may 
change  as  a  bifurcation  parameter  is  adjusted.  However,  throughout  the  range  of 
feedback  gains  from  8-14.5  mV/mW,  the  control  delay  that  produced  a  minimum 
in  the  error  signal  was  a  constant  value  of  19.315  ns  to  within  the  2.5  ps  resolution 
of  the  integration.  The  effects  of  varying  t>  will  be  investigated  in  more  detail  in 
the  next  chapter. 

To  locate  the  domain  of  control,  I  generate  many  time  series  like  those  in  Fig. 
4.9  using  a  range  of  values  of  7  and  7C  both  inside  and  outside  the  domain  of  control. 
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From  each  time  series,  I  calculate  the  root  mean  square  value  of  c(t)  as  defined  in 
Eq.  4.13.  Figure  4.10  shows  how  the  average  error  signal  changes  as  the  control 
gain  is  varied  at  three  different  values  of  7.  The  dotted  line  gives  {c(t))rms  as  a 
function  of  the  control  gain  with  7  =  8.52  mV/mW.  Near  7C  =  0  mV/mW,  the 
control  error  is  high  indicating  that  the  system  is  outside  the  domain  of  control. 
Around  7C  =  1.0  mV/mW,  (c(t))rms  abruptly  falls  two  orders  of  magnitude  to  a 
value  close  to  zero.  The  fact  that  it  does  not  go  all  the  way  to  zero  (within  machine 
precision)  is  due  to  the  fact  that  there  is  some  error  in  7>,  as  mentioned  earlier. 
Near  7C  =  6.0  mV/mW,  the  {c(t))rms  again  changes  abruptly,  this  time  increasing 
rapidly  by  two  orders  of  magnitude.  A  similar  behavior  is  displayed  by  the  system 
with  7  =  13.2  mV/mW.  The  long  dashed  line  in  Fig.  4.10  shows  {c(t))rms  in  this 
case.  At  two  points,  (c(t))rms  jumps  abruptly.  Between  these  points  the  error  signal 
is  very  small.  The  dotted  line  in  Fig.  4.10  shows  quite  different  behavior  in  the  case 
where  7  =  14.64  mV/mW.  Although  a  minimum  occurs,  the  decrease  is  less  than 
an  order  of  magnitude  and  no  abrupt  jump  occurs. 

I  interpret  those  regions  where  the  gain  is  extremely  low  as  being  within  the 
domain  of  control.  In  order  to  define  the  boundaries  of  the  estimated  domain,  a 
cutoff  value  of  {c(t))rms  must  be  chosen.  Since  the  error  signal  falls  so  rapidly  when 
control  is  achieved,  the  resulting  domain  estimate  should  not  depend  sensitively  on 
the  cutoff  chosen.  Based  on  the  appearance  of  time  series  from  borderline  values  of 
7c  I  set  the  cutoff  at  1CF2  mW.  The  resulting  estimated  domain  of  control  is  shown 
in  Fig.  4.11. 

The  estimated  domain  is  qualitatively  quite  similar  to  the  domain  estimated 
experimentally  (see  Fig.  4.8).  The  shape  is  triangular  showing  the  linear  decrease 
in  the  range  of  effective  control  gains  as  the  feedback  gain  is  increased.  Control  fails 
completely  near  7  =  14.5  mV/mW.  As  in  the  experiment,  this  is  the  beginning  of 
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Figure  4.11:  Domain  of  control  of  modified  TDAS  applied  to  the  model  of  the 
active  interferometer. 

the  region  where  the  dynamics  become  very  broad  band. 

4.4  Summary 

In  this  chapter,  I  have  proposed  a  modification  of  the  TDAS  control  method  suitable 
for  fast,  time-delay  chaotic  systems.  In  this  new  method,  the  control  loop  latency  is 
intentionally  set  to  equal  the  time  delay.  I  have  demonstrated  this  method  exper¬ 
imentally  by  applying  it  to  the  active  interferometer  with  bandpass  feedback  and 
external  modulation.  I  developed  a  model  of  the  system  with  control  that  produced 
behavior  in  good  agreement  with  the  experiment.  Two  time  delays  appear  in  the 
control  term  in  the  modified  TDAS  method.  In  the  next  chapter,  I  will  examine 
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the  effect  on  control  of  variation  of  these  delays. 
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Chapter  5 

Influence  of  Delay  M  ismatch  in  the 
M  odified  TDAS  Control  M  ethod 


The  modified  TDAS  control  scheme  introduced  in  the  preceeding  chapter  for  con¬ 
trolling  time-delay  systems  requires  adjusting  two  time  delays  in  the  control  loop. 
First,  the  control  perturbation  consists  of  the  difference  between  the  interferometer 
output  at  two  times  seperated  by  a  time  interval  which,  in  this  chapter,  I  will  refer 
to  as  the  control  time  delay.  Ideally,  the  control  delay  is  equal  to  the  period  of  the 
orbit  to  be  controlled,  7>.  The  second  delay  time  in  the  controller  is  the  control 
loop  latency,  which  is  set  to  equal  the  time  delay  td  of  the  feedback  loop  of  the 
uncontrolled  system.  The  modified  TDAS  control  method  was  motivated  by  the  ob¬ 
servation  that  many  time  delay  systems  have  a  strong  peak  in  the  auto  correlation 
function  around  td-  Since  this  peak  has  a  finite  width  and  is  not  centered  exactly 
on  td,  it  is  interesting  to  ask  how  varying  the  latency  around  Td  affects  control. 
Similarly,  how  does  the  varying  the  control  time  delay  away  from  the  true  orbit  of 
the  period  affect  control? 

Studies  of  the  original  TDAS  control  method  found  that  variation  of  these  delays 
produced  quite  distinct  effects  on  control.  Effective  TDAS  control  strongly  depends 
on  the  careful  matching  of  the  control  time  delay  with  the  period  of  the  unstable 
orbit  [49,  113].  For  example,  in  an  experimental  study  of  controlling  a  diode  res¬ 
onator  using  TDAS,  a  dramatic  increase  in  the  error  signal  was  observed  when  the 
control  time  delay  differed  by  a  few  hundred  picoseconds  from  the  orbit  period  of 
100  ns.  On  the  other  hand,  the  error  signal  remained  small  while  the  control  loop 
latency  was  adjusted  by  several  nanoseconds. 
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The  theoretical  explanation  for  these  differing  effects  is  as  follows.  With  the 
control  time  delay  exactly  equal  to  the  period  of  the  desired  orbit,  the  error  signal 
is  identically  zero  when  the  system  is  on  the  orbit.  Thus  the  unstable  orbit  of  the 
system  without  control  is  also  an  orbit  of  the  system  with  control.  Any  mismatch 
between  the  control  time  delay  and  the  orbit  period  results  in  a  non-zero  error  signal 
when  the  system  is  on  the  desired  orbit.  Therefore,  the  desired  orbit  of  the  system 
without  control  is  not  a  periodic  orbit  of  the  system  with  control  applied.  However, 
the  controlled  system  will  follow  a  periodic  orbit  that  differs  only  slightly  from  the 
desired  orbit  of  the  uncontrolled  system  if  the  mismatch  is  small  enough  [114].  This 
so-called  fictitious  orbit  has  a  shape,  period,  and  stability  that  tend  towards  those 
of  the  uncontrolled  orbit  as  the  mismatch  tends  toward  zero.  Since  the  period  of  the 
fictitious  orbit  differs  from  the  control  time  delay,  the  error  signal  tends  toward  zero 
as  the  mismatch  is  decreased.  Thus,  the  error  signal  is  expected  to  have  a  minimum 
when  the  control  delay  is  exactly  7>  and  to  increase  smoothly  as  the  mismatch  is 
increased. 

Variation  of  the  control  loop  latency  produces  a  different  response  from  the 
controller.  Latency  is  ideally  zero  in  TDAS  control,  but  experiments  have  shown 
successful  control  with  non- zero  latency  [49] .  Positive  control  loop  latency  does  not 
prevent  the  error  signal  from  going  to  zero  on  the  desired  orbit.  Therefore,  no  matter 
how  much  latency  is  present,  the  orbit  of  the  uncontrolled  system  is  always  an  orbit 
of  the  controlled  system  (assuming  the  control  delay  is  properly  set).  However,  the 
latency  does  affect  the  stability  of  the  orbit  [65,  66,  115].  Thus,  the  error  signal  is 
expected  to  be  zero  as  long  as  the  latency  is  within  some  stability  domain  but  to 
grow  abruptly  when  the  a  stability  boundary  is  crossed. 

The  modifications  to  TDAS  introduced  for  controlling  time-delay  systems  do 
not  alter  the  basic  role  played  by  the  control  delay  and  the  latency.  Therefore, 
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qualitatively  similar  results  are  expected  when  these  delays  are  varied  around  their 
nominal  values.  As  in  TDAS,  the  desired  orbit  of  the  uncontrolled  system  is  only 
an  orbit  of  the  system  with  control  if  the  control  time  delay  exactly  equals  7>.  The 
error  signal  should  have  a  unique  minimum  when  the  control  time  delay  equals  Tp. 
Also,  as  in  TDAS,  control  should  be  stable  in  a  finite  range  of  latency  around  the 
nominal  value  of  rL  =  Tp. 

To  test  these  expectations,  I  investigate  the  influence  of  variation  control  loop 
latency  and  control  delay  from  their  nominal  values  on  the  modified  TDAS  control 
method  as  applied  to  the  active  interferometer  with  delayed  bandpass  feedback.  In 
Sec.  5.1,  I  describe  simulations  using  the  mathematical  model  developed  in  the 
preceeding  chapters.  In  Sec.  5.2,  I  present  experimental  results  and  a  comparison 
with  the  model.  Both  the  experimental  and  simulated  results  are  consistent  with 
the  conclusions  drawn  above  regarding  the  role  of  the  delays  in  the  control  method. 
In  Sec.  5.3,  I  discuss  the  implications  of  these  results  on  the  practicality  of  the 
modified  TDAS  method. 

5.1  Simulations 

The  active  interferometer  with  delayed  bandpass  feedback  and  modified  TDAS  con¬ 
trol  is  described  by  the  equations 

tl V  =  -V  +  yP^  {1  +  &i  sin  [c at  (P^  -  P0)]}  +  (5.1) 

P  =  (P  —  P0)  +  k  V  +  flmVm  cos  ,  (5.2) 

th  L  -I 

c(t)  =  {Pm  {1  -  b2  sin  [a  (P^  -  P0)]}  -  P^  +Tp  {1  -  b2  sin  [a  (Pm  +Tp  -  P0)]}}  , 

(5.3) 

derived  in  Chapters  3  and  4.  The  quantity  c(t)  is  referred  to  here  as  the  error 
signal.  For  all  simulations  descibed  in  this  chapter,  I  set  the  feedback  loop  gain 
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(a)  (b) 


F  igure  5.1:  (a)  Experimental  and  (b)  theoretical  domains  of  control  as  determined 
in  the  previous  chapter.  The  black  circle  shows  the  point  in  each  domain  at  which 
the  data  in  this  chapter  is  collected. 

7  =  12.28  so  that  the  uncontrolled  system  displays  chaotic  oscillations  according 
to  Sec.  3.5.  The  control  gain  is  set  to  yc  =  4.65,  placing  the  system  within  both 
the  experimental  and  theoretical  domains  of  control  as  shown  in  Fig.  4.8.  All  other 
parameters  are  as  given  in  Ch.  4.  In  particular,  the  feedback  delay  td  =  19.1  ns. 

I  first  consider  variations  of  the  control  time  delay  and  fix  the  control  loop 
latency  at  77  =  19.1  ns.  I  numerically  integrate  Eqs.  5.1  and  5.2  using  nearly  the 
same  procedure  as  described  in  Sec.  4.3.  However,  I  use  a  time  step  ten  times 
smaller.  As  explained  in  Sec.  4.3,  the  integration  time  step  is  chosen  so  that  an 
integer  number  of  steps  equals  the  feedback  delay  time,  7y>  This  means  that,  in 
general,  the  period  of  the  orbit  7>  is  not  equal  to  an  integer  number  of  integration 
steps.  However,  as  the  integration  step  size  is  reduced,  an  integer  number  of  steps 
can  better  approximate  the  true  value  of  t>.  Using  a  time  steps  of  2.5  ps,  I  can 
approximate  the  true  orbit  period  to  within  less  than  0.02%  (assuming  the  true 
period  is  within  1  ns  of  Td). 

At  each  of  several  values  of  the  control  time  delay,  I  quantify  the  control  error 
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by  calculating  the  root  mean  square  average  of  c(t).  To  obtain  each  average  value,  I 


Controller  Delay  (ns) 


Figure  5.2:  Average  error  signal  in  the  modified  TDAS  controller  as  control  loop 
latency  is  varied  in  the  mathematical  model.  The  feedback  gain  7  =  12.28  and  the 
control  gain  7C  =  4.65. 

first  integrate  through  a  transient  of  duration  10,000td  and  then  average  the  value 
of  c(t)  over  the  next  1007y>  Figure  5.2  shows  (c(t))  on  a  logarithmic  scale  as  a 
function  of  control  time  delay.  As  expected,  a  very  sharp  minimum  is  observed. 
The  control  time  delay  is  equal  to  19.315  ns  at  the  minimum.  Adjusting  the  control 
time  delay  just  7.5  ps  above  or  below  this  value  results  in  an  increase  of  (c(f))  by  two 
orders  of  magnitude.  Adjusting  the  control  delay  by  an  additional  100  ps  above  or 
below  results  in  another  order  of  magnitude  increase  in  the  error  signal.  Finally,  at 
roughly  1  ns  to  either  side  of  the  19.315  ns,  ( c(t ))  is  four  orders  of  magnitude  above 
the  minimum  value.  The  sharp  minimum  and  smooth  increase  in  the  error  signal 
as  the  control  delay  is  moved  away  from  19.315  ns  is  consistent  with  the  theoretical 
picture  of  a  gradual  distortion  of  the  orbit  of  the  uncontrolled  system. 

Next,  I  investigate  the  influence  of  the  control  loop  latency  on  the  effectiveness 
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of  the  controller.  I  keep  the  value  of  the  control  time  delay  fixed  at  19.315  ns  and 
quantify  the  effectiveness  of  the  controller  by  determining  (c(t))  while  I  vary  the 
latency.  The  result  is  shown  in  Fig.  5.3.  Over  a  range  of  values  from  0  to  40  ns,  the 
average  value  of  c{t)  is  uniformly  high  except  for  a  deep  notch  strecthing  from  19.0 
ns  to  19.45  ns  where  the  error  drops  abruptly  by  four  orders  of  magnitude  to  the 
minimum  value  observed  in  the  previous  paragraph.  This  finite  region  where  the 
error  is  very  small  and  the  abrupt  transitions  at  the  edges  of  the  region  where  the 
error  signal  increases  dramatically  are  consistent  with  theoretical  description  given 
above  based  on  results  from  the  original  TDAS  method. 

It  is  notable  that  control  is  not  effective  in  the  model  when  the  latency  ap¬ 
proaches  zero.  A  theoretical  study  has  predicted  there  may  be  cases  where  ETDAS 
control  is  only  effective  in  some  non-zero  range  of  latency  but  no  example  has  ap¬ 
peared  in  the  literature  [115].  Also  notable  is  the  fact  that  control  fails  when  the 
latency  is  equal  to  2td .  Apparently,  at  least  at  this  control  gain  value,  the  mem¬ 
ory  of  the  system  decays  too  rapidly  for  successful  control  with  a  latency  equal  to 
multiples  of  td- 

5.2  Experimental  Observations 

In  this  section  I  describe  the  experimentally  observed  effects  of  varying  the  delays  in 
the  modified  TDAS  control  method  when  applied  to  the  active  interferometer  with 
delayed  bandpass  feedback.  The  experimental  setup  is  identical  to  that  described 
in  Sec.  4.2.  I  set  the  gain  7  =  12.28  so  that  the  uncontrolled  system  displays 
chaotic  oscillations.  The  control  gain  7C  =  4.65,  placing  the  system  well  within 
the  experimentally  determined  domain  of  control  shown  in  Fig.  4.8.  Following  the 
order  of  presentation  of  the  preceeding  section,  I  will  first  consider  variations  of  the 
control  delay  time  followed  by  variations  of  the  control  loop  latency. 
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18.6  18.8  19.0  19.2  19.4  19.6  19.8  20.0 
Control  Loop  Latency  (ns) 

Figure  5.3:  Average  error  signal  in  the  modified  TDAS  controller  as  control  loop 
latency  is  varied  in  the  mathematical  model.  The  feedback  gain  7  =  12.28  and  the 
control  gain  7C  =  4.65.  A  wide  range  of  latency  from  zero  to  over  2td  is  shown  in 
(a).  The  region  around  the  minimum  is  shown  in  more  detail  in  (b). 
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Figure  5.4:  The  control  delay  is  determined  by  the  difference  in  propagation  time 
through  the  two  highlighted  sections.  The  direct  path  is  shown  as  a  heavy  dashed 
line.  The  delayed  path  is  shown  as  a  heavy  solid  line. 

The  control  delay  is  determined  physically  by  the  difference  in  propagation  time 
through  the  two  paths  connecting  the  splitter  labelled  C  in  Fig.  4.5  and  the  hy¬ 
brid  coupler  labelled  E.  These  two  paths  are  shown  schematically  in  Fig.  5.4. 
Experimentally,  this  path  length  is  varied  by  adjusting  the  length  of  the  delay  line 
located  just  after  the  splitter.  Coarse  adjustments  are  made  by  adding  or  subtract¬ 
ing  fixed  lengths  of  RG-58  coaxial  cable.  A  20  cm  length  of  cable  adds  approxi¬ 
mately  1  ns  to  the  control  loop  latency.  Fine  adjustments  are  made  using  a  constant 
impedance,  variable  length  delay  line  with  a  tunable  range  of  310  ps  (General  Radio 
Co.  874-LK10L). 

To  quantify  the  effect  of  varying  the  control  delay,  I  record  the  root  mean  square 
average  of  the  error  signal.  A  high-speed  digital  oscilloscope  records  2,500  samples 
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Figure  5.5:  Average  error  signal  in  the  modified  TDAS  controller  as  the  control 
time  delay  is  varied.  The  circles  show  experimental  values.  The  feedback  gain 
7  =  12.28  and  the  control  gain  7C  =  4.65.  The  solid  line  shows  the  simulation 
reproduced  from  Fig.  5.2. 

of  the  error  signal  at  a  rate  of  1  GS/s  and  then  calculates  the  root  mean  square. 
In  Fig.  5.5,  the  circles  show  the  resulting  values  of  (c(t)).  The  solid  line  is  the 
simulated  result  from  Fig.  5.2  reproduced  on  a  linear  scale.  A  minimum  appears  at 
19.33  ns  where  the  error  signal  drops  to  the  noise  level  as  defined  in  Sec.  4.2.3.  The 
error  signal  rises  by  a  factor  of  two  within  0.5  ns  on  either  side  of  the  minimum. 
The  width  of  the  minimum  is  comparable  to  that  observed  when  the  control  loop 
latency  is  varied.  Although  no  sharp  minimum  is  apparent  as  in  the  simulations, 
the  gross  shape  of  the  experimentally  observed  minimum  agrees  with  the  simulated 
minimum. 

I  next  investigate  the  effect  of  varying  the  control  loop  latency.  In  the  active 
interferometer  with  bandpass  feedback,  the  control  loop  latency  is  the  total  propaga¬ 
tion  time  of  the  signal  as  it  travels  through  the  closed  loop  formed  by  the  controller, 
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Figure  5.6:  Schematic  diagram  of  active  interferometer  with  moclifed  TDAS  con¬ 
trol.  The  control  loop  latency  is  the  propagation  time  around  the  path  shown  here 
by  the  heavy  black  line. 

laser  and  interferometer.  This  path  is  shown  schematically  in  Fig.  5.6  by  the  heavy 
black  lines.  The  path  length  is  varied  by  adjusting  the  length  of  the  delay  line 
located  just  after  the  detector.  As  in  the  previous  experiment,  a  variable  delay  is 
included  for  fine  adjustments.  I  set  the  control  time  delay  to  19.33  ns,  the  location 
of  the  observed  minimum  in  Fig.  5.5. 

The  result  is  shown  in  Fig.  5.7.  The  error  signal  displays  a  minimum  at  19.1  ns 
reaching  the  estimated  noise  level  of  ~  4  mV  (see  Sec.  4.2.3).  It  roughly  triples  in 
size  1  ns  above  or  below  the  minimum.  The  transition  from  the  region  where  the 
error  signal  is  small  to  the  region  where  it  is  large  is  somewhat  smoother  than  in  the 
simulated  result.  However,  the  minimum  is  clearly  much  wider  than  was  observed 
when  the  control  time  delay  was  varied  as  predicted  above. 
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Figure  5.7:  Average  error  signal  in  the  modified  TDAS  controller  as  the  control 
loop  latency  is  varied.  The  circles  show  experimental  values.  The  feedback  gain 
7  =  12.28  and  the  control  gain  7C  =  4.65.  The  solid  line  shows  the  simulation 
reproduced  from  Fig.  5.3. 

5.3  Discussion 

In  this  chapter,  I  examined  the  role  of  the  control  time  delay  and  control  loop 
latency  in  the  modified  TDAS  control  method.  The  effect  on  the  control  error  signal 
when  each  of  these  quantities  was  varied  around  its  nominal  value  was  observed  in 
simulations  and  in  an  experiment.  The  success  of  control  was  found  to  depend  more 
sensitively  on  the  control  time  delay  than  on  the  latency.  The  accuracy  required  in 
setting  these  delays  is  a  central  consideration  in  determining  the  practicality  of  the 
method. 

My  observations  suggest  that  the  latency  needs  to  be  known  only  to  within 
several  percent  of  its  true  value  to  achieve  control  in  many  situations.  Any  value 
within  this  finite  range  results  in  successful  control  with  a  very  small  error  signal. 
Therefore,  it  does  not  appear  to  be  a  significant  barrier  to  the  practical  application 
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of  the  method.  The  fact  that  the  range  of  latencies  over  which  control  is  successful  is 
located  around  td  rather  than  around  zero  in  the  original  TDAS  method  constitutes 
the  primary  advantage  of  the  modified  technique.  Roughly  speaking,  a  controller 
based  on  the  original  TDAS  method  applied  to  a  fast  time-delay  system  will  require 
faster  components  than  those  used  in  the  implementation  of  the  system  itself.  On 
the  other  hand,  a  controller  based  on  the  modified  TDAS  method  can  be  constructed 
with  similar  components  to  those  used  in  the  system  itself. 

The  situation  is  more  complicated  with  regard  to  the  control  time  delay.  The 
error  signal  is  very  sensitive  to  this  parameter  and  goes  to  zero  only  at  when  the 
control  time  delay  is  precisely  equal  to  the  period  of  the  uncontrolled  orbit.  In  sys¬ 
tems  such  as  the  active  interferometer  examined  here,  the  presence  of  noise  obscures 
the  very  sharp  minimum  in  the  error  signal.  The  usefulness  of  the  technique  for  a 
noisy  system  then  rests  on  whether  the  broader  minimum  in  the  error  signal  that 
is  observable  is  deep  enough  for  the  application  at  hand.  On  the  other  hand,  in 
systems  that  are  subject  to  much  lower  noise  levels,  the  sharp  minimum  makes  the 
control  method  a  useful  tool  for  accurately  determining  the  period  of  the  unstable 
orbits  of  a  chaotic  system  similar  to  schemes  devised  for  automating  the  process  of 
locating  the  correct  period  in  the  original  TDAS  method  [116,  117]. 
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Chapter  6 
Conclusion 


The  ability  to  control  fast  chaotic  systems  is  an  important  tool  for  building  a  funda¬ 
mental  understanding  of  the  dynamics  of  these  systems  and  for  exploiting  them  in 
technological  applications.  This  dissertation  documents  my  investigation  of  control¬ 
ling  a  fast  chaotic  time-delay  system.  To  conclude  this  report,  I  review  the  major 
results  reported  in  each  chapter  and  suggest  further  experiments  and  applications. 

6.1  Summary 

The  first  two  chapters  of  this  thesis  provide  introductory  and  background  material. 
In  Chapter  1,  I  delineate  the  place  of  my  research  in  the  larger  fields  of  nonlinear 
dynamics  and  physics.  This  general  introduction  is  followed  by  a  chapter-by-chapter 
preview  of  the  thesis.  In  Chapter  2,  I  provide  a  more  specific  context  for  my  work. 
I  define  the  concept  of  control  loop  latency  and  illustrate  it  with  a  simple  example 
of  proportional  control  of  a  linear  system.  I  show  that  the  domain  of  control  for  this 
simple  system  is  dramatically  reduced  as  the  control  loop  latency  is  increased.  I 
then  review  the  development  of  techniques  for  controlling  chaos  in  general  and  fast 
chaos  in  particular. 

In  Chapter  3,  I  introduce  a  new  fast  chaotic  opto-electronic  device  consisting 
of  a  semiconductor  laser,  a  Mach-Zehnder  interferometer,  and  a  bandpass  filtered 
delayed  electronic  feedback  loop.  The  device  offers  a  high  degree  of  design  flexi¬ 
bility  at  a  cost  much  lower  than  other  known  sources  of  fast  optical  chaos.  The 
nonlinearity  in  the  system  is  determined  by  the  path  difference  through  the  arms  of 
the  interferometer.  The  time  scale  of  the  chaotic  dynamics  is  determined  by  a  the 
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delay  time  and  bandwidth  of  the  feedback  loop.  Thus,  both  of  these  key  character¬ 
istics  of  the  system  are  easily  accessible  and  manipulable.  For  example,  by  simply 
adjusting  the  feedback  delay  time  in  my  experimental  system,  the  fundmental  fre¬ 
quency  of  oscillation  can  be  set  anywhere  in  range  30-103  MHz.  In  principle,  the 
system  could  easily  be  modified  to  oscillate  at  1-10  GHz.  This  flexibilty  makes  the 
system  suitable  for  laboratory  studies  of  dynamics  at  slow  speeds.  I  characterize 
the  behavior  of  the  system  as  the  feedback  loop  gain  is  varied  both  when  the  system 
is  running  autonomously  and  when  it  is  subject  to  external  modulation.  I  develop 
a  mathematical  model  of  the  system  that  reproduces  the  experimentally  observed 
behavior.  In  later  chapters  of  this  thesis,  the  system  serves  as  a  testbed  for  a  new 
method  for  controlling  fast  time-delay. 

In  Chapter  4,  I  propose  a  new  method  for  controlling  chaos  in  a  fast  time-delay 
system.  The  method  is  a  modification  of  the  TDAS  control  technique  where  the 
control  loop  latency  is  precisely  set  equal  to  td,  the  feedback  time  delay  in  the 
uncontrolled  system.  I  demonstrate  this  technique  by  applying  it  to  the  active  in¬ 
terferometer  with  bandpass  feedback  described  in  the  previous  chapter.  The  domain 
of  control  is  estimated  experimentally  and  numerically.  The  experimental  and  the¬ 
oretical  results  show  good  qualitative  agreement.  The  domain  is  found  to  have  a 
rather  simple,  triangular  shape  that  pinches  off  at  a  high  feedback  loop  gain.  Prior 
to  my  work,  the  fastest  controlled  chaotic  system  reported  in  the  literature  had  a 
characteristic  oscillation  frequency  of  19  MHz  with  a  control  loop  latency  of  4.4  ns 
[63].  With  the  modified  TDAS  method,  control  is  achieved  when  the  fundamental 
frequency  of  oscillation  is  51.8  MHz  depite  the  much  larger  latency  of  19.1  ns. 

In  Chapter  5,  I  investigate  the  role  of  the  control  time  delay  and  control  loop 
latency  in  the  modified  TDAS  control  method.  Results  from  studies  of  the  time 
delays  in  the  original  TDAS  method  are  reviewed  and  their  application  to  the  mod- 
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ified  TDAS  method  is  discussed.  I  examine,  through  experiments  and  simulations, 
the  effect  of  varying  the  control  time  delay  and  latency  around  their  nominal  values. 
The  observed  bahvior  is  consistent  with  expectations  based  on  the  studies  of  the 
original  TDAS  method.  The  control  time  delay  must  be  set  exactly  equal  to  the 
period  of  the  desired  UPO  if  the  control  perturbations  are  to  vanish  when  the  sys¬ 
tem  is  actually  on  the  orbit.  In  contrast,  the  control  loop  latency  may  be  adjusted 
within  a  finite  range  without  affecting  the  success  of  control. 

6.2  Future  Directions 

The  two  main  contributions  of  this  thesis  to  the  field  of  nonlinear  dynamics  are 
the  introduction  of  a  flexible  new  source  of  fast  optical  chaos  and  a  new  method 
for  controlling  fast  time-delay  chaotic  systems.  Both  of  these  innovations  open  the 
door  to  further  investigations  and  applications. 

Much  work  remains  to  be  done  to  fully  characterize  the  dynamics  of  the  active 
interferometer  with  delayed  bandpass  feedback.  In  Chapter  3,  I  observed  a  wide 
range  of  behavior  as  the  feedback  gain  was  varied.  A  more  complete  picture  would 
be  provided  by  a  careful  study  and  classification  of  the  bifurcations  that  occur 
between  the  first  periodic  oscillation  and  broad  band  chaos.  It  is  also  important 
to  understand  the  complexity  of  the  dynamics  both  as  a  function  of  the  feedback 
gain  and  the  time  delay.  Only  a  couple  of  studies  of  the  dimension  of  systems 
with  bandpass  delayed  feedback  have  been  reported.  In  both  cases,  the  estimated 
dimension  was  surprisingly  high  relative  to  similar  systems  with  lowpass  delayed 
feedback  [44,  87]. 

The  active  interferometer  with  delayed  bandpass  feedback  has  many  features 
that  may  make  it  a  good  candidate  for  uses  in  chaos-based  communications  ap¬ 
plications.  First,  the  use  of  a  semiconductor  laser  facilitates  its  integration  with 


117 


modern  optical  communication  technologies  such  as  optical  fiber  transmission,  op¬ 
tical  switching,  etc.  The  visible  wavelength  of  the  diode  laser  (660  nm)  used  here  was 
chosen  for  experimental  convenience  over  the  invisible  infra-red  wavelengths  (1.3- 1.5 
/jrn)  commonly  used  in  fiber  optic  communications  systems.  However,  redesign  at 
communications  wavelengths  would  be  straight  forward. 

Second,  chaos  synchronization  may  be  easier  to  achieve  with  this  system.  In 
general,  high-quality  synchronization  requires  nearly  identical  chaotic  systems  [118]. 
Important  parameters  of  the  active  interferometer  with  delayed  bandpass  feedback, 
such  as  the  nonlinearity  and  the  delay  time  are  accessible  and  easily  fine  tuned.  The 
main  potential  obstacle  to  high-quality  synchronization  in  this  system  is  the  phase 
noise  of  the  diode  laser  [119]. 

Third,  the  use  of  external  electrical  modulation  offers  a  means  for  encoding  mes¬ 
sages  by  techniques  successfully  demonstrated  using  much  more  expensive  chaotic 
sources  such  as  the  erbium-doped  fiber  laser.  The  basic  idea  is  to  use  the  message  as 
a  driving  signal  for  the  chaotic  system.  The  message  gets  wrapped  up  in  the  chaotic 
dynamics  in  a  complicated  way  making  extraction  practically  impossible  without  a 
synchronized  reciever  [27]. 

Another  interesting  application  of  the  active  interferometer  may  be  the  rapid 
generation  of  random  numbers  [120,  121].  A  wide  variety  of  computing  applications 
rely  heavily  on  random  numbers.  For  example,  random  numbers  are  used  as  encryp¬ 
tion  keys  to  protect  data  from  unauthorized  decryption.  The  ever  increasing  speed 
of  modern  processors  produces  a  need  for  physcial  systems  that  can  generate  the 
truly  random  numbers  at  fast  and  faster  rates.  A  study  of  the  statistical  properties 
of  the  oscillations  observed  in  the  active  interferometer  is  necessary  to  determine 
whether  the  system  passes  the  standard  tests  for  randomness. 

The  modified  TDAS  control  method  provides  a  new  tool  for  the  investigation 
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and  exploitation  of  fast  chaotic  systems  like  the  active  interferometer  with  bandpass 
delayed  feedback.  In  this  thesis,  a  proof-of-concept  experiment  demonstrates  the 
successful  application  of  the  technique.  However,  much  more  can  be  done  to  under¬ 
stand  how  and  when  it  works.  Theoretical  methods  for  predicting  the  domain  of 
control  with  more  accuracy  than  simple  numerical  integration  already  exist  for  the 
original  TDAS  method  [64],  An  extension  of  these  methods  to  the  modified  TDAS 
method  is  desirable.  Also,  the  modified  TDAS  method  could  be  experimentally 
tested  with  a  more  well-known  chaotic  system  such  as  an  electronic  circuit  imple¬ 
mentation  of  the  Mackey- Glass  system  [122].  The  dependence  of  the  dimension  of 
this  system  on  the  time  delay  is  well  characterized  and  thus  the  effectiveness  of  this 
method  in  controlling  high-dimensional  chaos  could  be  more  precisely  determined. 
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